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ABSTRACT 

We have performed an unbiased deep near-infrared survey toward the Aquila molecular 
cloud with a sky coverage of ~1 deg^. We identified 45 molecular hydrogen emission-line ob- 
jects(MHOs), of which only 11 were previously known. Using the Spitzer archival data we also 
identified 802 young stellar objects (YSOs) in this region. Based on the morphology and the 
location of MHOs and YSO candidates, we associate 43 MHOs with 40 YSO candidates. The 
distribution of jet length shows an exponential decrease in the number of outflows with increasing 
length and the molecular hydrogen outflows seem to be oriented randomly. Moreover, there is 
no obvious correlation between jet lengths, jet opening angles, or jet H 2 1 — 0 ^(l) luminosities 
and spectral indices of the possible driving sources in this region. We also suggest that molecular 
hydrogen outflows in the Aquila molecular cloud are rather weak sources of turbulence, unlikely 
to generate the observed velocity dispersion in the region of survey. 

Subject headings: stars: formation - stars: winds, outflows - ISM: jets and outflows - infrared: ISM - 
shock waves 


1. Introduction 

Mass outflow plays an essential role in the pro¬ 
cess of star formation (Shu et al. 1987; Arce et al. 
2007; Bally et al. 2007). It is believed that mass 
outflow is an important way to transfer the ex¬ 
cess angular momentum from the dense molecular 
cores to the ambient interstellar medium (Shang 
et al. 2007). Outflows have been observed in dif¬ 
ferent wavelengths: Herbig-Haro (HH) objects are 
the optical manifestation of shock-ionized mass 
outflows, which trace either material ejected from 
the protostars directly or the shocked interstel¬ 
lar medium (Reipurth & Bally 2001; Wang et al. 
2004; Walawender et al. 2005; Wang et al. 2005; 
Wang & Henning 2006, 2009). CO outflows de¬ 
tected in millimeter wavelength probe the swept- 
up or entrained medium along the edges of the jets 
(Bachiller 1996; Wu et al. 2004). 


In near-infrared bands, molecular hydrogen 
emission lines, in particular the v =1-0 S(l) tran¬ 
sition at 2.12 /im, are powerful tracers of shock- 
excitation. Davis et al. (2010) has defined molecu¬ 
lar hydrogen emission-line objects (MHOs) as the 
near-infrared manifestation of mass outflow and 
presented a comprehensive catalog of over 1400 
MHOs. Wide-held deep near-infrared survey us¬ 
ing a combination of narrow-band H 2 Alter and the 
corresponding Alter for continuum emission such 
as the Kg Alter has become a very efAcient method 
to detect and identify MHOs (Stanke et al. 2002; 
Khanzadyan et al. 2004; Davis et al. 2008; Froe- 
brich et al. 2011). MHOs are also good tracers of 
young stellar objects (YSOs), in particular those 
deeply embedded in the molecular cloud cores. 
Once the relationship between the MHOs and 
their driving sources has been established, we can 
study the distribution of YSOs, the interaction of 
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YSOs with the ambient interstellar medium, the 
star formation efficiency, and the evolution of the 
entire star-forming region by analyzing the sta¬ 
tistical characteristics of the molecular hydrogen 
outflows (Stanke et al. 2002; Davis et ah 2008, 
2009; Zhang et al. 2013). 

The Aquila Rift is located along the Galactic 
plane and the large Aquila Rift cloud complex 
stretches from 20° to 40° in longitude and -1° to 
10° in latitude, as revealed by CO and H I obser¬ 
vations (Dame et al. 2001; Prato et al. 2008). The 
molecular mass of the Aquila Rift estimated from 
CO observation is about 1.1-2.7x 10^ Mq (Dame & 
Thaddeus 1985; Dame et al. 1987; Straizys et al. 
2003). Prato et al. (2008) reviewed the Aquila 
Rift and compiled a list of known YSOs that in¬ 
cludes only 9 YSOs from literature in the Aquila 
Rift. The apparent lack of a large number of YSOs 
in Aquila is especially surprising given the age of 
young stars detected in this cloud (Prato et al. 
2003; Rice et al. 2006), and the abundance of raw 
material for star-formation in this region. How¬ 
ever, in the recent past, with the availability of 
data from the Spitzer, WISE, and Herschel space 
telescopes, a number of new YSOs have been de¬ 
tected in this region. 

The Aquila Rift has now become a hot-spot 
for the study of star formation. Gutermuth et al. 
(2008b) discovered an embedded cluster of YSOs 
in the Serpens-Aquila Rift using the Spitzer IRAC 
imaging data and identified 54 Class I and 37 
Class II sources in the cluster. Similarly, as 
one of the targets of the Herschel Gould Belt 
key program (Andre et al. 2010), the Aquila re¬ 
gion has also been surveyed with the PACS and 
SPIRE (Konyves et al. 2010; Bontemps et al. 2010; 
Men’shchikov et al. 2010). Bontemps et al. (2010) 
used “the Aquila Rift complex” to describe the 
molecular cloud that corresponds to a large ex¬ 
tinction feature in the extinction map derived us¬ 
ing the 2MASS catalog. The Aquila Rift com¬ 
plex mainly harbors two known sites of star for¬ 
mation (see Fig. 1 in Bontemps et al. 2010): Ser¬ 
pens South is the western young embedded cluster 
(Gutermuth et al. 2008b) and W40 is the east¬ 
ern cluster associated with an H II region (Smith 
et al. 1985; Vallee 1987). “The Aquila molecu¬ 
lar cloud” referred to in this paper is the same as 
the “Aquila Rift complex” described by Bontemps 
et al. (2010). The distance estimates to the Aquila 


Rift complex vary from ~200pc to up to 900 pc 
(Radhakrishnan et al. 1972; Vallee 1987; Straizys 
et al. 2003; Gutermuth et al. 2008b; Rodney & 
Reipurth 2008; Shaping et al. 2012). Bontemps 
et al. (2010) compared the different estimates of 
distance and finally adopted the value of 260 pc. 
In this paper, we follow their suggestion and also 
adopt the value of 260 pc as the distance of the 
Aquila molecular cloud. 

Hundreds of protostars have been discovered 
in the Aquila molecular cloud (Bontemps et al. 
2010). Nakamura et al. (2011) identified 15 
blueshifted and 10 redshifted outflow components 
using the GO (J = 3-2) mapping observations in 
the Serpens South cloud. Connelley et al. (2007) 
identified several H 2 features around IRAS 18264- 
0143 in the Aquila region. Teixeira et al. (2012) 
identified 14 MHOs and associated them with 10 
YSOs in the Serpens South cloud. These works 
suggest that star-format ion in the Aquila molec¬ 
ular cloud is still in its youth which makes this 
cloud an ideal laboratory to study various phases 
of stellar birth. 

In this paper we present our wide-held near- 
infrared survey toward the Aquila molecular cloud 
with a coverage of ~1 square degree. Our goal is 
to detect the jets and outflows from young stars 
in the Aquila molecular cloud. Compared with 
the optical observations, our near-infrared survey 
enables us to detect more outflow features due 
to the relatively lower extinction in near-infrared 
bands. Combining our near-infrared observations 
of this region with the corresponding archival data 
from the Spitzer and Herschel, we were able to as¬ 
sociate H 2 emission features with their possible 
driving sources. We then studied the statistical 
characteristics of H 2 outhows in the Aquila molec¬ 
ular cloud. This paper is organized as follows 
: In Section 2 we describe our near-infrared ob¬ 
servations and the complementary archival data. 
We then present our results in Section 3, includ¬ 
ing the scheme employed to identify outhows and 
YSOs. In Section 4, we discuss the characteris¬ 
tics of molecular hydrogen outhows in the Aquila 
molecular cloud, including a statistical analysis of 
outhow parameters such as jet lengths, jet open¬ 
ing angles, jet orientations and the H 2 (1-0) S(l) 
luminosity of jets. We also calculate the momen¬ 
tum injected by molecular hydrogen outhows in 
this cloud and critically evaluate its contribution 
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to turbulence that appears to support the Aquila 
molecular cloud against self-gravity. Finally we 
summarize our findings in Section 5. 

2. Observations and data reduction 
2.1. Near-infrared imaging 

The observations were conducted in queue- 
scheduled observing (QSO) mode between 26th 
and 29th July, 2012 with WIRCam (Puget et al. 
2004) equipped on Canada-France-Hawaii Tele¬ 
scope (CFHT), covering in total an area of ~1 
deg^. WIRCam is a near-infrared mosaic im¬ 
ager with four 2048x2048 CCDs, yielding a field 
of view of ~21'x21' with a plate scale of Oi'SOG 
pixel” We observed 10 fields toward the Aquila 
molecular cloud in J, iJ, Kg, and H 2 bands. The 
narrowband H 2 filter has a central wavelength of 
2.122/xm and a bandwidth of 0.032/xm^. Each 
field was imaged with a four point dithering pat¬ 
tern. The number of sub-exposures per dithering 
pattern position for J, H, Kg, and H 2 band is 
2, 3, 2, and 3, respectively, while the individual 
exposure time is 30 s in J band, 10.3 s in H band, 
15 s in Kg band, and 100 s in H 2 band. Thus the 
total integration time of J, H, Kg, and H 2 band 
for each field is 240s, 123.6s, 120s, and 1200s, re¬ 
spectively. Figure 1 shows the coverage of our 
CFHT/WIRCam observations with red boxes. 

Individual images are primarily processed by 
the CFHT ‘I’iwi pipeline (the IDL Interpreter of 
the WIRCam Images, version 2.1.1), which in¬ 
cludes dark subtraction, flat-fielding, non-linearity 
correction, cross-talk removal, and sky subtrac¬ 
tion^. Then, the data are handled with the SIM¬ 
PLE IDL package^, which is an IDL based data re¬ 
duction package for optical and near-IR blank-field 
imaging observations. We use SIMPLE-WIRCAM 
(SIMPLE Imaging and Mosaicking PipeLinE for 
WIRCAM) to remove the distortion and do the 
absolute astrometry on a single frame by compar¬ 
ing the image with the 2MASS reference catalog 


^The information of the broadband filters can be found 
at http://www.cfht.hawaii.edu/Instruments/Filters/ 
wircam.htinl 

^The details of the pipeline can be found at 
http://www.cfht.hawaii.edu/Instruments/Imaging/ 
WIRCam/IiwiVersionlDoc.html 
^http://www3.asiaa.sinica.edu.tw/-whwcuig/idl/SIMPLE/ 
index.htm 


(Skrutskie et al. 2006). To estimate the accuracy 
of our astrometry, we use Sextractor (Bertin & 
Arnouts 1996) to do source detection in each fil¬ 
ter of each field and compare the coordinates of 
detected sources with the WCS entries of 2MASS 
catalog. We find that the root-mean-square (rms) 
values of the difference of coordinates between the 
sources on the final mosaic images and 2MASS 
catalog for all fields are below 0.3". Based on the 
accurate astrometry, the dithered individual expo¬ 
sures are finally combined into the stacked images 
with the software SWARP (Bertin et al. 2002), 
which is a program that resamples and co-adds 
together PITS images using any arbitrary astro- 
metric projection defined in the WCS standard^. 

Point source detection and aperture photome¬ 
try® is performed on the mosaic images via the 
IDL routine find and aper with an aperture ra¬ 
dius of 1.5" and concentric sky annuli of inner and 
outer radii of 3.0" and 4.5" respectively. Pho¬ 
tometric zero points and color terms were cal¬ 
culated by comparison of the instrumental mag¬ 
nitudes of relatively isolated, unsaturated bright 
sources with the counterparts in the 2MASS Point 
Source Catalog (Skrutskie et al. 2006). By com¬ 
paring ~23000 sources detected in both the CEHT 
and the 2MASS observations, we find that the 
photometric reliability in all bands is ~0.08-0.16 
mag, depending on the source brightness. 

2.2. Spitzer data 

As part of the Gould Belt Legacy program 
(PID: 30574), the Spitzer Space Telescope ob¬ 
servations toward the Serpens-Aquila rift were 
conducted in May and September, 2007 with 


^http://www.astromatic.net/software/swarp 
®As our survey region is close to the Galactic plane, the 
more reliable PSF fitting method for source detection and 
photometry is still under construction. By comparing the 
results from aperture photometry and PSF fitting for one 
field that is closest to the Galactic plane, we found that 
the mean value and root mean square value of difference 
of photometry between these two methods are <0.02 mag 
and <0.07 mag individually for all bands, but that PSF 
fitting method can detect larger number of faint sources. 
In this paper, photometry for the point sources in our near- 
infrared images is only used to select YSO candidates as the 
supplement of Spitzer data in the YSO selection scheme 
suggested by Gutermuth et al. (2009). Thus the aperture 
photometry is enough to match our science objective. The 
PSF fitting results and follow-up analysis will be published 
in our subsequent paper (Sun et al. 2015, in preparation). 
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the IRAC and MIPS cameras (Fazio et al. 2004; 
Rieke et al. 2004). IRAC images at 3.6, 4.5, 5.8, 
and 8.0 /rm were made in High Dynamic Range 
mode with integration times of 0.4 and 10.4 s. 
We download the IRAC standard basic calibrated 
data (BCD) products provided by the Spitzer Sci¬ 
ence Center from their standard data processing 
pipeline version S18.I8. Final mosaics are built at 
the native instrument resolution of 1.2" pixel“^ for 
each Astronomical Observation Request (AOR), 
using Mopeofi (version 18.5, Makovoz & Khan 
2005), which is a package for reducing and an¬ 
alyzing imaging data. Thus for each AOR, we 
obtain two IRAC mosaics, one built from the long 
exposures of 10.4 s and the other built from the 
short exposures of 0.4 s. MIPS images at 24, 70, 
and 160 ^m were obtained at the fast scan rate of 
17" s“^. We also download the MIPS post-BCD 
data products provided by the Spitzer Science 
Center from their pipeline version S18.12. Note 
that we only use the MIPS 24 pja images with 
the resolution of 2.45" pixel”^ in this paper. The 
Spitzer observations toward the Serpens-Aquila 
rift cover a very large sky region. To match with 
our near-infrared imaging data, we finally restrict 
the IRAC and MIPS data to the region that is 
marked with the black dashed box in Fig. 1. 
Figure 2 (left panel) shows the three-color im¬ 
age that is constructed with IRAC 3.6/rm (blue), 
4.5y^m (green), and 8.0/rm (red) images. 

Point source detection and aperture photome¬ 
try is performed on the final IRAC and MIPS mo¬ 
saics via the photometry and visualization tool, 
PhotVis (version 1.10; Gutermuth et al. 2004, 
2008a). Aperture photometry is performed on the 
IRAC images with an aperture radius of 2.4" and 
background flux is estimated with an concentric 
sky annuli of inner and outer radii of 2.4", and 
7.2" respectively. Aperture and inner and outer 
sky annulus radii are selected to be 7.6", 7.6", and 
17.8", respectively, for MIPS 24/rm images. Pho¬ 
tometric zero points for IRAC aperture photome¬ 
try are derived from the calibrations presented in 
Reach et al. (2005), including standard aperture 
corrections for the radii adopted. The photomet¬ 
ric zero point of MIPS 24 /rm band is the suggested 
value from Gutermuth et al. (2008a). The 90% 


®http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/mopex/ 


completeness limits are derived by adding succes- 
sively dimmer sets of Spitzer PRFs (Point Re¬ 
sponse Functions^), extracting their fluxes using 
the same procedure described above. The detec¬ 
tion completeness limits are defined as the mag¬ 
nitudes at which 90% of the synthetic stars are 
recovered. Our 90% completeness limits are 14.7, 
15.1, 13.4, 13.0, and 8.2 mag at 3.6, 4.5, 5.8, 8.0, 
and 24 /tm. 

2.3. Multi-band photometric catalog 

To obtain the final multi-band photometric cat¬ 
alog, we performed the bandmerge in four stages. 
First, we crossmatch our CFHT photometric cat¬ 
alog with the 2MASS point source catalog (Skrut- 
skie et al. 2006) using a tolerance of 2". The pho¬ 
tometry of bright sources from our CFHT pho¬ 
tometric catalog is contaminated by saturation. 
Thus for the point sources brighter than J = 14.5, 
or 77 = 13.5, or Kg = 13 mag we use 2MASS 
photometric results. Second, we combine the long 
exposure photometry and short exposure photom¬ 
etry for each IRAC band and then the four IRAC 
band source lists are matched with a tolerance of 
2". Third, the MIPS 24/rm catalog is integrated 
with a wider tolerance of 4" due to the low reso¬ 
lution (6") of the MIPS 24/rm band. Finally, we 
crossmatch the Spitzer photometric catalog and 
the near-infrared photometric catalog with a tol¬ 
erance of 2". Note that the final catalog uses the 
mean IRAC and MIPSl positions for all entries. 

2.4. Herschel extinction map 

The Herschel archival data used in this paper 
are part of the Herschel Gould Belt guaranteed 
time key programmes for the study of star forma¬ 
tion with the PACS (Poglitsch et al. 2010) and 
SPIRE (Griffin et al. 2010) instruments (Andre 
et al. 2010) and have been published in Andre 
et al. (2010), Konyves et al. (2010), and Bontemps 
et al. (2010). The details of the observations and 
data reduction can be found in Konyves et al. 
(2010). We use the final calibrated images® at 160, 


^see http://irsa.ipac.caltech.edu/data/SPITZER/ 

docs/irac/calibrationfiles/psfprf/ and http: 
//irsa.ipac.caltech.edu/data/SPITZER/docs/mips/ 
calibrationfiles/prfs/ 

®http://www.herschel.fr/cea/gouldbelt/en/Phocea/Vie_ 
des_labos/Ast/ast_visu.php?id_ast=66 
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250, 350, and 500 /im with the angular resolutions 
of ~12", ~18", ~25", and ~37", respectively. 

All Herschel maps are smoothed to the beam 
size of the 500 /rm map, through convolving the 
maps with the convolution kernels® supplied by 
Gordon et al. (2008). The column density was 
determined from a pixel-to-pixel modified black 
body ht to four longer wavebands of PACS and 
SPIRE from 160 /rm to 500 /rm, assuming the dust 
opacity law of Beckwith et al. (1990) and dust 
emissivity index of 2 (Hildebrand 1983), follow¬ 
ing the same procedure described by Konyves 
et al. (2010). Note that we do not use the PACS 
70 /xm map to construct the column density map 
because 70/xm data may not be entirely tracing 
the cold dust, as suggested by Hill et al. (2011). 
The column density map is calibrated with the 
2MASS extinction map from Dobashi (2011). The 
final Herschel extinction map is obtained using 
the relation Abxj = 10^^ x Ay (Bohlin et al. 
1978; Konyves et al. 2010). Figure 2 (right panel) 
shows the obtained Herschel extinction map of the 
Aquila molecular cloud, overlaying the contours 
with levels of Ay = 4, 10, and 30 mag. 

3. Results 

3.1. The detected MHOs in Aquila 

For each field, we have obtained the JI 2 narrow- 
band image and the continuum Kg band image. 
We use Sextractor (Bertin & Arnouts 1996) to 
detect point sources on H 2 images and Kg im¬ 
ages. Then the unsaturated bright point sources 
are used to adjust the flux level of stars in the H 2 
image and Kg image and a continuum-subtracted 
H 2 image {H 2 — Kg) is obtained for each field. 

We did the visual inspection on all H 2 — Kg im¬ 
ages to search for H 2 line emission features.These 
H 2 line emission features are treated as MHO can¬ 
didates. In order to avoid the inclusion of instru¬ 
mental artifacts, each MHO candidate is examined 
in the corresponding H 2 , J, H, and Kg images. 
Any potential artifact has been removed from the 
MHO candidate list. 

Emission nebulae, such as planetary nebulae 
(PNe), supernova remnants (SNRs) and H II re¬ 
gions may contaminate our identification. PNe 


®http: //dirty. as. arizona.edu/~kgordon/niips/conv_ 
psf s/conv_psf s.html 


usually exhibits symmetrical morphology and can 
be distinguished from MHOs based on their mor¬ 
phological difference. We found no PNe in our 
MHO candidate list. Moreover, the known H II 
region, W40, is located in our survey region of the 
Aquila molecular cloud. We examined MHO can¬ 
didates near W40 and two (^2%) potential con¬ 
taminants have been removed from our MHO can¬ 
didate list. The remaining extended H 2 emission 
features are identified as MHOs. 

Finally, we have identified 45 MHOs that con¬ 
sist of 108 discrete H 2 emission-line features. All 
are situated within the Aquila molecular cloud 
(see Fig. 3). Of these 45 MHOs, 11 are previously 
known objects while 34 are newly discovered. Ta¬ 
ble Al in the appendix lists the position, surface 
area, radius, and integrated H 2 line flux for each 
sub-feature within each MHO. The MHO numbers 
assigned to each object have been entered into the 
MHO catalog^® (Davis et al. 2010). A description 
of each MHO is also given in the appendix. 

3.2. Photometry of MHO features 

Areal photometry is performed to measure the 
fluxes of MHO features. We define manually a 
polygonal aperture based on the morphology and 
surface brightness distribution of each MHO fea¬ 
ture. The principle for the aperture definition is 
to ensure that no stars are in the aperture and the 
aperture contains as little background area outside 
the bounds of each H 2 emission feature as possible. 
Note that there are circular holes in some polyg¬ 
onal apertures in order to avoid stars. Figs. Al- 
A42 in the appendix show these apertures with 
blue polygons and the holes are also marked with 
red circles. For each MHO feature, a same poly¬ 
gon aperture toward the nearby sky region that is 
emission free is used to estimate the local sky back¬ 
ground. The polygon apertures around the H 2 
features and the nearby sky apertures are applied 
to the MHO features on the continuum-subtracted 
H 2 images to measure the fluxes of the H 2 1-0 S(l) 
emission line. 

The flux calibration is obtained via our near- 
infrared photometric catalog (2MASS-I-CFHT, see 
section 2.3). We use Sextractor (Bertin & Arnouts 
1996) to measure the fluxes of the point sources 
on the H 2 images. The zero-point flux is derived 


^®http://astro.kent.ac.uk/-df/MHCat/ 


5 



through the comparison of the integrated counts 
with the cataloged flux of the point sources. We 
then applied the zero-point flux to the integrated 
counts of the MHO features with the aim to con¬ 
vert the counts of MHO features to the flux unit. 
Uncertainties in the areal photometry are esti¬ 
mated from the variation of the local sky back¬ 
ground level in the continuum-subtracted H 2 im¬ 
ages. 

The obtained fluxes and uncertainties of the 
MHO features are shown in table A1 in the ap¬ 
pendix. In our H 2 images, outflow features with a 
surface brightness of ~1.3xl0“^® W m“^ arcsec“^ 
are detected at ~5 (t above the surrounding back¬ 
ground. The median of the fluxes of MHO fea¬ 
tures is 1.2 X 10“^^ W m“^ and the minimum of 
the fluxes we detected is 6.6 x 10“^^ Wm“^. 

3.3. Identification of YSOs in Aqnila 

YSOs usually show excessive infrared emission 
that can be used to distinguish them from field 
stars and classify them into different evolutionary 
classes. Lada & Wilking (1984) and Lada (1987) 
developed an empirical classification scheme that 
firstly codified a tripartite class system based on 
the slope of the SED to indicate the different evo¬ 
lutionary stages of YSOs, which has become a 
‘standard’ classification system of YSOs. Some 
authors use the infrared spectral index to iden¬ 
tify YSO candidates in the star-forming regions 
(Mallick et al. 2013; Kim et al. 2015). YSOs can 
also be classified using their mid-IR clours (i.e., 
color-magnitude diagrams). Gutermuth et al. 
(2008a, 2009) established a mid-IR color-based 
method to identify YSOs, which has been widely 
applied to the YSO identification in the star¬ 
forming regions. We use this scheme in the present 
work as it can efficiently mitigate the effects of 
contamination from field stars and extragalactic 
sources. We calculate the slope of SEDs and esti¬ 
mate the evolutionary stages of the YSOs detected 
in this survey using the classification scheme sug¬ 
gested by Lada (1987). 

YSO identification is based on our multi-band 
photometric catalog (see section 2.3) and the 
source classification scheme presented in Guter¬ 
muth et al. (2008a, 2009). The details about this 
multiphase source classification scheme can be 
found in the appendix of Gutermuth et al. (2009). 
Here we just summarize our process (also see the 


description in Rapson et al. (2014)). 

There are three phases in the YSO selection 
scheme suggested by Gutermuth et al. (2009). 
Phase 1 is applied to the sources that have de¬ 
tections in all four IRAG bands with photometric 
uncertainties cr <0.2 mag. After removing the con¬ 
taminants such as star-forming galaxies, broad¬ 
line AGNs, unresolved knots of shock emission, 
and sources that have PAH-contaminated aper¬ 
tures, the sources with [4.5]-[5.8] > 0.7 mag and 

[3.6] -[4.5] >0.7 mag are classified as Glass I candi¬ 
dates while the sources with [4.5]-[8.0] > 0.5 mag, 

[3.6] -[5.8] > 0.35 mag, [3.6]-[4.5] > 0.15 mag, and 

[3.6]-[5.8] < 3.5x ([4.5]-[8.0]-0.5)-|-0.5 mag are con¬ 
sidered as Class H candidates, accounting for pho¬ 
tometric uncertainty. The remaining sources are 
classified as Class Hl/field sources. 

Phase 2 is applied to the sources that lack de¬ 
tections at either 5.8 or 8.0 /im, but have high 
quality (cr <0.1 mag) near-infrared detections in 
J, H, and Ks bands. In order to distinguish 
the sources with IR-excess from those that are 
simply reddened by dust along the line of sight, 
we deredden the photometry of sources based on 
the extinction law presented in Wang & Jiang 
(2014) and Elaherty et al. (2007). The dereddened 
ATs—[3.6] and [3.6]-[4.5] colors are used to iden¬ 
tify the sources with infrared excess at 3.6 /rm and 
4.5 /im accounting for photometric uncertainty. 
Sources with no IR excess at 3.6 /xm and 4.5 /xm are 
presumed as Class Hl/field sources. 

Phase 3 is applied to the sources that have 
detections in MIPS 24 /xm band with the pho¬ 
tometric uncertainties cr <0.2 mag. The Class 
Hl/field sources that were classified in previous 
two phases are re-examined and sources with col¬ 
ors of [5.8]-[24] > 2.5 mag or [4.5]-[24] > 2.5 mag 
and [3.6] < 14 are classified as transition disk can¬ 
didates. Eor the sources that lack detection in 
some IRAG bands, but are very bright at 24/xm, 
i.e., [24] < 7 mag and [IRACX]-[24] > 4.5 mag, 
they are classified as deeply embedded Class I pro¬ 
tostars, where [IRACX] is the photometric magni¬ 
tude for the longest wavelength IRAG detection. 
The AGN candidates and shock emission domi¬ 
nated sources that were classified as contaminants 
in phase 1 are also re-examined and identified as 
Class I candidates if they have both bright MIPS 
24/xm photometry ([24] < 7 mag) and convinc¬ 
ingly red IRAC/MIPS colors of [3.6]-[5.8] > 0.5 
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mag and [4.5]-[24] > 4.5 mag and [8.0]-[24] > 4 
mag. Finally, the Class I source that were classi¬ 
fied in all three phases are re-analyzed to ensure 
their nature. Sources with [5.8]-[24] > 4 mag or 
[4.5]-[24] > 4 mag are classified as true Class I 
sources while other sources that do not match the 
above requirements are identified as highly red¬ 
dened Class II sources. Figure 4 (left panel) shows 
the color-color diagrams that display the results 
of applying the above classification criteria to our 
multi-band photometric catalog. 

Finally we have identified 802 YSO candidates 
using the above criteria (Gutermuth et al. 2009), 
including 151 Class I sources, 629 Class II sources, 
and 22 transition disks. Figure 2 (left panel) shows 
the spatial distribution of these YSO candidates. 

The source spectral index, the slope of the 
source’s spectral energy distribution (SED) that 
is defined as 

dlog(XS(X)) 

(X = - 

dlog{X) 

where S{X) is the source’s flux density at wave¬ 
length A, can be used to determine the evolution¬ 
ary state of a source. We calculate the spectral 
indices of all 802 YSOs through fitting their ob¬ 
served SEDs from 2 /im to 24 /rm. Using the YSO 
classification scheme suggested by Lada (1987), we 
reclassify these 802 YSOs into 193 Class I sources, 
601 Class II sources, and 8 Class III/MS sources. 

The Aquila molecular cloud is highly obscured. 
The mean visual extinction in this region esti¬ 
mated from our Herschel extinction map is ^5 
magnitude. Thus for the YSO candidates obtained 
above, it is desirable to correct their flux densities 
for extinction. However, due to the lack of spec¬ 
tral types of YSO candidates, we can only roughly 
determine their extinction. 

For the YSO candidates with detections in J, 
H, and bands, we adopt the method suggested 
by Fang et al. (2013), in which the extinction is 
obtained by employing the H — Kg versus J — H 
color-color diagram. The detailed description of 
this method can be found in Fang et al. (2013). 
Here we just summarize some only a few aspects 
of this scheme. The location of each YSO in the 
H — Kg versus J — H color-color diagram de¬ 
pends on both its intrinsic colors and its extinc¬ 
tion. Figure 4 shows the H — Kg versus J — H 
color-color diagram of the YSO candidates in the 


Aquila molecular cloud. Given the different ori¬ 
gins of intrinsic colors of YSOs, we divide the 
color-color diagram into three sub-regions. Differ¬ 
ent methods are used to obtain the intrinsic color 
of the YSOs in different regions: in region 1, the 
intrinsic color [ J — H]o values of YSOs are simply 
assumed to be 0.6, which is the typical value of 
a K5-type dwarf star; in region 2, the [J — H]q 
value of a YSO candidate is obtained from the in¬ 
tersection between the reddening vector and the 
locus of main-sequence stars from Bessell & Brett 
(1988); in region 3, the intrinsic [J — H] color is 
derived from where the reddening vector and the 
CTTS locus (Meyer et al. 1997) intersects. Then 
the extinction values of individual YSOs are es¬ 
timated from their observed color [J — H]obs and 
their intrinsic color [J — H]o, using the equation 
Aif = {[J — H]obs — [J — H]o)/{^j/^K — AH/AK), 
where the extinction ratios were adopted from 
Wang & Jiang (2014). 

For the YSO candidates outside these three re¬ 
gions or without detection in J, H, or Kg bands, 
their extinction is estimated from our Herschel ex¬ 
tinction map. We firstly obtain the mean values 
of nine beam-size (~37" per beam) pixels around 
the YSO candidates and then use the half of these 
mean values as the extinction values of individ¬ 
ual YSO candidates, assuming that all the YSO 
candidates are located in the middle of the Aquila 
molecular cloud along the line of sight. Here we 
must stress that it is very difficult to accurately 
derive the extinction towards YSOs even if their 
spectral types are known because of the contribu¬ 
tion from the disk and accretion flow. Thus the 
specific value of extinction towards some YSO can¬ 
didates obtained using the above methods could 
have a large uncertainty. 

For colour excesses in the near-infrared we 
adopted the extinction law suggested by Wang 
& Jiang (2014), while adopting the one suggested 
by Flaherty et al. (2007) for colour excesses in the 
mid-infrared. The median value of visual extinc¬ 
tion toward the YSO candidates is 11 magnitudes. 
Based on the de-reddened flux densities of 802 
YSO candidates, we calculate their de-reddened 
spectral indices and reclassify them according to 
their de-reddened spectral-index, a. Finally, we 
obtain 100 Class I sources, 551 Class H sources, 
and 151 Class HI sources. Figure 5 shows the dis¬ 
tributions of apparent and de-reddened spectral 
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indices of 802 YSO candidates. 

3.4. Driving sources of MHOs in Aquila 

The likely driving sources of MHOs were iden¬ 
tified on the basis of the morphologies and loca¬ 
tions of individual MHOs relative to the candi¬ 
date YSOs using which, we associated 43 MHOs 
with 40 YSO candidates that led us to identify 
40 molecular hydrogen outflows. Table 1 lists the 
information of these 40 molecular hydrogen out¬ 
flows. Figure 3 shows the spatial distribution of 
the molecular hydrogen outflows with blue lines 
that connect the MHOs and their likely driving 
sources. Note that there are two MHOs (MHO 
3284-3285) whose driving sources can not be iden¬ 
tified due to the lack of YSOs around them. The 
details about the driving source identification can 
be found in the appendix. 

Of 40 outflow sources, 30 (75%) are protostars 
(a >0) when we use the apparent spectral indices 
to classify the YSOs. After correcting fluxes of 
outflow sources for extinction, 27 (~68%) of 40 
outflow sources belong to protostars. Thus ~70% 
H 2 outflows in the Aquila molecular cloud are 
driven by protostars. This fraction is lower than 
the value of >80% in Orion A (Davis et al. 2009) 
and higher than the value of <50% in Ophiuchus 
(Zhang et al. 2013). The mean value of the appar¬ 
ent spectral indices of outflow sources for Aquila 
is ~0.6, which is lower than the value of 0.86 in 
Orion A (Davis et al. 2009) and higher than the 
value of -0.16 in Ophiuchus (Zhang et al. 2013). 
We note that the detection limit for H 2 features 
(~1.3xl0“^® W m“^ arcsec“^) of our survey to¬ 
ward the Aquila molecular cloud is better than 
that toward Orion A (7x10“^® W m“^ arcsec“^, 
Davis et al. 2009) and similar to that toward Ophi¬ 
uchus (1x10“^® W m“^ arcsec“^, Zhang et al. 
2013). The difference in the mean value of the 
apparent spectral indices is likely due to various 
selection effects discussed by Zhang et al. (2013). 
With better sensitivity, we detect more low power 
jets from Class H sources in the Aquila molecu¬ 
lar cloud than in Orion A. As the Aquila molec¬ 
ular cloud is located at a distance of ~260 pc, 
farther than Ophiuchus (120 pc, Lombardi et al. 
2008; Loinard et al. 2008), less outflows with small 
spatial extent will be detected in Aquila than in 
Ophiuchus when the surface brightness detection 
limits are similar. 


We can also estimate that about 16% of proto¬ 
stars in our survey region drive H 2 outflows before 
de-reddening. After de-reddening for YSOs, the 
fraction of protostars that drive H 2 outflows turns 
out to be 27%; this fraction is respectively, 30% for 
the Orion A (Davis et al. 2009), and 20% for the 
Ophiuchus molecular clouds (Zhang et al. 2013). 
We also find that, of our 40 outflow sources, 18 are 
associated with infrared nebulae. Figure 6 shows 
the As—band images of these 18 infrared nebulae 
Note that several outflow sources are invisible in 
near-infrared. 

3.5. Interesting individual objects 

3.5.1. A disk-jet system 

Figure 7 shows the region in the vicinity of the 
MHO 3251bl and its driving source ID^17, which 
constitute a disk-jet system. The left panel is the 
H 2 image and the right panel is the Kg image. The 
solid line in the left panel shows the collimated 
jet, whose bright head has been identified as MHO 
3251bl (see Fig. A7 in Appendix). It seems that 
there are several faint knots to the northwest of 
MHO 3251bl along the solid line, but these struc¬ 
tures are too faint to be identified with any degree 
of certainty. Another obvious structure in Fig. 7 
is the bipolar infrared nebula. A dark lane is also 
visible between the two bright lobes and the ori¬ 
entation of the dark lane is roughly perpendicular 
to the jet. The white plus in the right panel marks 
the position of ID^17, which is the mean position 
of IRAC and MIPS 24 /im weighted with flux in 
individual band. We also mark the mean position 
of IRAC for ID#17 using the blue cross and the 
MIPS 24 /rm position using the red cross. We note 
that all positions of ID#I7 are roughly located in 
the dark lane, indicating a protostar surrounded 
by an edge-on disk. The apparent size of the dark 
lane at the base of the northern lobe of the neb¬ 
ula is ~4", corresponding to ~1000 AU for an as¬ 
sumed distance of 260 pc. The physical size of the 
bipolar nebula is ~I800 x 1300 AU, which is ~2 
times larger than those for the detected edge-on 
systems in Taurus (Padgett et al. 1999). 

3.5.2. An H 2 outflow from a Class 0 source 

Figure 8 shows the region of MHO 3260 (also 
see Fig. AI2) which is an H 2 outflow from ID#2 
that is identified as a Class 0 source by Maury 



et al. (2011). The top-left panel is a three-color im¬ 
age constructed with CFHT J (blue), H 2 (green), 
and Kg (red) images and MHO 3260 exhibits as 
green features in the panel. The detailed de¬ 
scription about MHO 3260 can be found in the 
appendix. The top-right panel shows the three- 
color image constructed with IRAC 3.6^m (blue), 
4.5/im (green), and 8.0/rm (red) images. The 
IRAC 4.5 /im band contains H 2 {i' = 0 — 0, 
S(9, 10, 11)) lines and CO (ii = 1 — 0) band 
heads (Reach et al. 2006) and has proved an ef¬ 
fective tool for exhibiting the mid-infrared emis¬ 
sion from outflows. Some green features that 
correspond to MHO 3260 can be seen in this 
panel. These extended green objects (EGOs) rep¬ 
resent the shock-excited emission at 4.5 micron, 
which are usually suggested as mid-infrared out¬ 
flows (Cyganowski et al. 2008; Chambers et al. 
2009; Zhang & Wang 2009). The bottom-left 
panel is the three-color image constructed with 
MIPS 24/im (blue), PACS 70/im (green), and 
160 /im (red) images while the bottom-right panel 
is the three-color image constructed with SPIRE 
250/im (blue), 350/im (green), and 500/im (red) 
images. We note that is invisible at 1.2- 

24 /im and visible in Herschel 70-500 /tm bands. 
Actually, Maury et al. (2011) also detected the 
counterpart of ID^2 in 1.2 mm dust continuum 
image and classified it as a Class 0 source. Thus 
MHO 3260 is an H 2 outflow driven by a Class 0 
source. 

4. Discussion 

4.1. Outflow statistics 

Based on the positions and photometry of MHO 
features, we calculate important outflow param¬ 
eters such as jet length (L), jet opening angle 
(0), jet position angle (PA), and jet H 2 I — 
0 S'(l) line luminosity according to the suggestions 
from Davis et al. (2009); loannidis & Froebrich 
(2012a,b); Zhang et al. (2013). 

A jet opening angle (9) is measured from a cone 
with the smallest vertex angle centered on the 
driving source that includes all of the detected H 2 
line emission associated with each MHO outflow 
in our continuum-subtracted H 2 image. Thus the 
opening angle can be obtained even for the outflow 
that consists of only one feature. A jet position 
angle (PA) is measured east of north as computed 


from the bisector of the opening angle. For the 
bipolar outflows, we select the angles smaller than 
180 degrees as the position angles. A jet length 
(L) is defined as the maximal distance from the 
pixels of MHO features to the source projected to 
the bisector of the opening angle. For the bipolar 
outflows, we sum up both lengths over two lobes. 
We have obtained the flux for each MHO feature in 
section 3.2. The flux of each molecular hydrogen 
outflow is obtained through summing up all fluxes 
over the MHO features in the outflow. Then the 
H 2 1 — 0 S'(l) luminosity of each outflow is calcu¬ 
lated using the relation of Lh 2 = Fh^ x dTrd^, as¬ 
suming the distance of 260 pc for the Aquila molec¬ 
ular cloud. Note that we did not consider the ef¬ 
fect of extinction. Table 1 lists these parameters 
and Fig. 9 shows the distribution of these outflow 
parameters. Note that these parameters are only 
calculated with MHOs (not considering the infor¬ 
mation provided by HH objects and CO outflows) 
and uncorrected for inclination to the line of sight. 

4 . 1 . 1 . Jet lengths 

We have obtained the projected lengths for 
all 40 molecular hydrogen outflows in the Aquila 
molecular cloud (see table 1). The distribution 
of jet lengths shows a steep decrease in the num¬ 
ber of outflows with increasing length in Fig. 9. 
In our sample, there are two outflows (5%) longer 
than Ipc, 30 outflows (75%) longer than 0.1 pc. 
It seems that the fraction of parsec scale outflows 
in our sample (5%) is similar to that of ~8-9% in 
Orion A (Stanke et al. 2002; Davis et al. 2009), but 
smaller than ^12% in the west of Perseus (Davis 
et al. 2008). The median length of our sample is 
0.24 pc, with the maximum of 1.65 pc. 

Although the distribution of jet lengths in our 
sample can be fitted with both, an exponential 
and a power-law, we find that the distribution is 
closer to an exponential behavior than a power-law 
function. Shown in Fig. 10 is a histogram with log¬ 
arithmic axes of the jet lengths; the solid line here 
represents a linear fit to this distribution. The 
dotted lines show the linear fittings for the distri¬ 
butions of jet lengths in different bin sizes. The 
slopes obtained from the different fittings are from 
-0.66 to -1.03, depending on the bin size. Here we 
adopt a mean value of slopes to describe the be¬ 
havior of number N of outflows with flow length 
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in the following way: 

N oc 

where the uncertainty of slope is the standard de¬ 
viation of all fitting slopes. This value of the expo¬ 
nent, —0.85, is consistent with the value obtained 
by loannidis & Froebrich (2012b). loannidis & 
Froebrich (2012a,b) detected 134 MHOs in the 
Galactic plane toward Serpens and Aquila using 
the UKIRT telescope and identihed 131 molecular 
hydrogen outflows. They fitted the distribution of 
jet lengths with the exponential function and ob¬ 
tained a slope of —0.75. The outflows detected by 
loannidis & Froebrich (2012a,b) are not associated 
with the Aquila molecular cloud due to the large 
distances of >2kpc of outflows in their sample. 
Their survey traces more luminous outflows and 
about 25% of outflows in their sample are parsec- 
scale flows. Thus it seems that outflows from low- 
mass stars and massive/intermediate stars may 
have the similar behavior in projected lengths. 

Figure 11 (left two panels) shows the relation 
between jet lengths and the spectral indices of 
driving sources. Note that outflow source ID^2 
is invisible at 2-24 /xm while the outflow source 
1D^19 is saturated in all IRAC bands and MIPS 
24 /xm band. To plot these sources on the figure, 
we assume the value of 3 for their spectral indices 
because they are associated with the protostars 
identified by Bontemps et al. (2010). Two longest 
outflows are both driven by Class I sources. Ac¬ 
tually, the mean length of outflows from Class 1 
sources is 0.34 pc while the mean length of out¬ 
flows from Class Il/Ill sources is 0.22 pc if using 
the apparent spectral indices to do the YSO clas¬ 
sification. After we correct the spectral indices for 
extinction, these two values become 0.36 pc and 
0.22 pc. It seems that Class I sources drive longer 
flows than Class II/III sources. However, from 
Fig. 11 we can not find any obvious correlation 
between jet lengths and spectral indices of driving 
sources (correlation coefficients < 0.3). Of course, 
we can not expect a simple linear correlation be¬ 
tween jet lengths and spectral indices of driving 
sources because the relation between them should 
be complex: for the very young sources, we should 
expect a short jet length as a jet should start from 
the length of zero; then jets probably expand, and 
later on get shorter again as force support fades 
(Stanke 2000; Zhang et al. 2013). However, we 


can neither see any non-linear relation between 
jet lengths and spectral indices of driving sources 
from Fig. 11. Similar results have been found in 
other star-forming regions such as Orion A (Davis 
et al. 2009), Perseus (Davis et al. 2008), and Ophi- 
uchus (Zhang et al. 2013). The molecular hydro¬ 
gen outflows in Orion A, Perseus, and Ophiuchus 
also show no correlation between jet lengths and 
outflow source spectral indices. 

4-. 1-2. Jet opening angles 

Figure 9 (top right) shows the distribution of jet 
opening angles of outflows detected in the Aquila 
molecular cloud. The median opening angle of our 
outflow sample is 31°, with a minimum of 1.3° and 
a maximum of 149°. 

Jet opening angle can be used to describe the 
collimation of an outflow. Some millimeter ob¬ 
servations indicate that protostars that are in an 
earlier stage of evolution have a higher likelihood 
of driving more collimated CO outflows (Lee et al. 
2002; Arce & Sargent 2006). However, this trend 
has not been found in the molecular hydrogen out¬ 
flow samples. Davis et al. (2008, 2009) and Zhang 
et al. (2013) investigated the relation between H 2 
jet opening angles and spectral indices of driving 
sources in Perseus, Orion A, and Ophiuchus, re¬ 
spectively, and they all found no obvious corre¬ 
lation between jet opening angles and the spec¬ 
tral indices of driving sources. Figure 11 (middle 
two panels) shows the relation between jet open¬ 
ing angles and spectral indices of driving sources 
in our outflow sample. The correlation coefficients 
between jet opening angles and apparent or de- 
reddened spectral indices of driving sources are 
both <0.3. However, the mean opening angle of 
outflows driven by Class I sources in our sample is 
48° while the mean opening angle of outflow driven 
by Class H/Hl sources is 22° if using apparent 
spectral indices as the YSO classification criteria. 
If we use the de-reddened spectral indices to do 
the YSO classification, the mean opening angles of 
outflows driven by Class I sources and Class H/Hl 
sources become 51° and 20°, respectively. Despite 
the lack of a linear correlation between jet open¬ 
ing angles and outflow source spectral indices, it 
seems that Class I sources in the Aquila molec¬ 
ular cloud drive H 2 outflows with larger opening 
angles. Does this mean that Class I sources have 
more likelihood to drive poor collimated H 2 out- 
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flows? We must note two facts: first, the H 2 out¬ 
flow with the largest opening angle is driven by 
a Class II source; second, the H 2 outflows driven 
by Class II sources in our sample may be very 
likely incomplete. Zhang et al. (2013) estimated 
that the fraction of Class II sources which drive 
H 2 outflows is ~15% in Ophiuchus, but this frac¬ 
tion in our sample is <2%. Thus we cannot es¬ 
timate the collimation of H 2 outflows driven by 
Class II source in the Aquila molecular cloud us¬ 
ing our incomplete sample. If we only consider 
the H 2 outflows driven by protostars (a >0) in 
the Aquila molecular cloud, their opening angles 
show no obvious correlation with the spectral in¬ 
dices of their driving sources (correlation coeffi¬ 
cients <0.3), being consistent with the studies in 
other star-forming regions (Davis et al. 2008, 2009; 
Zhang et al. 2013). 

Davis et al. (2009) discussed the reasons for the 
lack of correlations between jet lengths or jet open¬ 
ing angles and spectral indices of driving sources. 
Firstly, shock-excited H 2 emission is not a good 
tracer of outflow parameters due to its short cool¬ 
ing time. H 2 emission is not a sensitive tracer for 
long jets, which results in the under-estimation of 
jet lengths. In addition, the wings of jet-driven 
bow shocks are often wider than the underlying 
jets and changes in flow direction due to preces¬ 
sion, which results in the over-estimation of jet 
opening angles. Secondly, the precise relationship 
between source spectral index and source age has 
not been established yet. 

4-. 1.3. Jet position angles 

Previous studies of Orion A (Stanke et al. 2002; 
Davis et al. 2009), Perseus (Davis et al. 2008), and 
Ophiuchus (Zhang et al. 2013) found that the ori¬ 
entation of outflows shows a homogeneous distri¬ 
bution with no significant trends. However, the 
study of DR21/W75 by Davis et al. (2007) found 
that the molecular hydrogen outflows, in particu¬ 
lar from massive cores, are preferentially orthog¬ 
onal to the molecular ridge, indicating a physi¬ 
cal connection between PAs of outflows and the 
large-scale cloud structure. A histogram of the 
position angles of molecular hydrogen outflows de¬ 
tected in the Aquila molecular cloud is shown in 
Fig. 9 (bottom left panel), using a bin size of 30°. 
Note that here we assume all outflows to be bipo¬ 
lar outflows and transfer their position angles to 


the range of [0°, 180°]. It seems that there is a 
peak at 60°-90°, but the value of the peak is only 
several counts higher than nearby columns. Actu¬ 
ally, a Kolmogorov-Smirnov (KS) test shows that 
there is a >93% probability that such a distribu¬ 
tion is drawn from a homogeneously distributed 
sample. Thus here we propose the distribution 
of jet position angles to be a homogeneous distri¬ 
bution, which means that outflows in the Aquila 
molecular cloud may orientate randomly. 

4 . 1 . 4 . Jet H 2 1-0 S(l) luminosities 

The H 2 1 — 0 ^(l) luminosities of our out¬ 
flows range from ~3x 10“® to ^2.6x 10“^ Lq with 
the median of 4.7x10“® L© (see table 1). Fig¬ 
ure 9 (bottom right panel) shows the distribution 
oi H 2 1 — 0 S'(l) luminosities of H 2 outflows in 
logarithm (in bins of 0.5). This distribution ex¬ 
hibits a decrease in the number of outflows with 
increasing H 2 1 — 0 S'(l) luminosities in the range 
of log{L 2 .i 2 /Lq) > —4. Note that there is also a 
decrease in the number of objects with very low 
A 2.12 H 2 outflows (log(L 2 .i 2 < -5)), which could 
be due to the incompleteness of low luminosity 
H 2 outflow sample. 

Stanke et al. (2002) and loannidis & Froebrich 
(2012b) investigated the outflow luminosities in 
Orion A and Galactic plane toward Serpens and 
Aquila, respectively. The 772 1 — 0 «5'(1) luminosi¬ 
ties of outflows in Orion A range from ^10“^ to 
^10-2 Lq while the luminosities of outflows in the 
Galactic plane toward Serpens and Aquila range 
from ^10“® to ~1O“^L0. They also And that 
the H 2 1 — 0 S'(l) luminosities of outflows exhibit 
a power-law behavior. The slopes of power-law 
functions obtained by them are —1.1 in Orion A 
and —1.9 in the Galactic plane toward Serpens and 
Aquila, respectively. For our sample, the distribu¬ 
tion of H 2 1-0 S'(l) luminosities is not a well- 
peaked histogram. The peak of this distribution 
changes with the histogram bin size. Gonsequently 
we did not statistically test the distribution of jet 
H 2 1 —0 S'(l) luminosities in the Aquila molecular 
cloud. 

Stanke (2000) investigated the evolutionary 
trends in the flow H 2 luminosities in Orion A and 
found that more evolved flow sources seem to drive 
lower luminosity outflows. Khanzadyan et al. 
(2012) also found a relatively poor correlation 
between outflow luminosity in H 2 1 — 0 S'(l) line 
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emission and spectral indices of driving sources 
based on 13 outflows in Braid Nebula. Figure 11 
(right two panels) shows the relation between out¬ 
flow H 2 1—0 S'(l) luminosities and spectral indices 
of driving sources in the Aquila molecular cloud. 
We find no obvious correlation (correlation coef¬ 
ficient < 0.2) between H 2 1 — 0 S')!) luminosities 
of our outflows and spectral indices of driving 
sources. The mean H 2 1 — 0 S(l) luminosity of 
outflows from Class I sources is SxlO^^L© while 
the mean H 2 1 — 0 S(l) luminosity of outflows 
from Class II/III sources is IxlO^^L© if using 
de-reddened spectral indices to do the YSO clas¬ 
sification. Thus although no clear correlation is 
found, it seems that there is still a general trend 
that Class I sources drive more luminous H 2 out¬ 
flows than Class II/III sources. Of course, this 
conclusion is not robust. At best, it is only a mod¬ 
est observation. However, even if this conclusion 
were true, we should like to point our readers to 
two caveats : (i) the contribution oi H 2 1 — 0 S') 1) 
line to the H 2 luminosity of outflows is only ~ 10 % 
(Caratti o Garatti et al. 2006) and, (ii) H 2 fluxes 
are affected by the environment-Class I sources 
are more likely to be in dense clouds where there 
is more H 2 gas that can be shock excited into 
emission while Class II sources are probably more 
widely distributed in regions where there isn’t 
very much H 2 gas. Therefore, any strong sug¬ 
gestion about the evolution of outflow luminosity 
merely encouraged by observations of the H 2 (I- 0 ) 
S(l) line-emission is likely to be too far-fetched. 

4.2. Momentum injection and cloud sup¬ 
port 

Feedback from mass outflows injects turbu¬ 
lent energy and momentum into the parent cloud 
that could possibly support the cloud against self¬ 
gravity. In this section we discuss the energy 
budget of the turbulence injected by the detected 
molecular hydrogen outflows in the Aquila molecu¬ 
lar cloud. To this end we calculate the momentum 
that is likely to be injected by this sample of out¬ 
flows and investigate if this injected momentum is 
sufficient to maintain the observed turbulence in 
the Aquila molecular cloud. 

Following the suggestions by Walawender et al. 
(2005); Davis et al. (2008); loannidis & Froe- 
brich ( 2012 b), we estimate the turbulent momen¬ 
tum in the Aquila molecular cloud using its mass 


(Me) and velocity dispersion {SVtur), as Ptur ~ 
MeSVtur- 

We use our Herschel extinction map obtained 
in section 2.4 to estimate the cloud mass of our 
survey region. Assuming the gas-to-dust ratio 
of N{H 2 )/Av =0.94x 10^^ cm“^ mag“^ (Bohlin 
et al. 1978) and a distance of 260 pc for the Aquila 
molecular cloud, the total mass in our survey re¬ 
gion is about 4200 Mq. However, the momentum 
injection from outflows is only likely in the rela¬ 
tively dense regions, i.e., those regions having col¬ 
umn density above the star-formation threshold 
for the cloud. Lada et al. (2010) found that the 
star formation rate (SFR) in molecular clouds is 
linearly proportional to the cloud mass above an 
extinction threshold of Ay ~ 7.3 mag, correspond¬ 
ing to a gas volume density threshold of n{H 2 ) « 
10^cm“^. Adopting this density as the threshold 
above which star-formation in the Aquila molec¬ 
ular cloud is likely, we estimate that the mass of 
potentially star-forming gas in this region is ~2700 
Mq. 

From Fig. 2, we can see that there are mainly 
two clouds in our survey region, one associated 
with the Serpens South cluster and another asso¬ 
ciated with the H II region, W40. We mark the 
location of W40 in Fig. 2 and Fig. 3 with circles 
whose scales are in proportion to the size of W40. 
The size of W40 (7') is estimated using radio con¬ 
tinuum observation data by Mallick et al. (2013). 
The distance estimation of W40 varies from ~300 
to ~900pc by different authors (Radhakrishnan 
et al. 1972; Vallee 1987; Rodney & Reipurth 2008; 
Shuping et al. 2012). Recent study by Shuping 
et al. (2012) suggested a distance of ~500pc for 
W40 based on the near-infrared spectral analysis 
of OB stars in the W40 H II region, which indicates 
that W40 may be not associated with Serpens and 
Serpens South cloud. We also note that, of our 40 
outflows, only one may be associated with W40. 
After excluding the mass of W40 (mass inside the 
circle marked in Figs. 2 and 3), the remain mass 
above star formation threshold is about 2300 M©. 

The turbulent velocity is difficult to measure in 
the molecular cloud. Here we follow the suggestion 
from Davis et al. (2008); loannidis & Froebrich 
( 2012 b) and use the widths of molecular lines to 
estimate the turbulent motions. The typical value 
of volume density {n{H 2 )) for the brightest NH 3 
emission in dark clouds is ~ 10^—10^cm“^ (Ho 
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& Townes 1983), indicating that NH 3 can be a 
tracer of dense molecular cloud with star forma¬ 
tion activity. Levshakov et al. (2013) presented 
the results of mapping observations in the NH 3 ( 1 , 
1) and ( 2 , 2 ) lines toward 49 sources in the Aquila 
region with the Effelsberg 100 m telescope. They 
detected NH 3 emission lines in 19 sources and 
the widths of lines are from ~ 0.2 to 1.5kms“^. 
We adopt a mean value of 0.7kms“^ to estimate 
turbulent velocity SVtur- Note that this value is 
also similar to the mean width of C ^®0 ( 2 - 1 ) lines 
( 0 . 6 kms“^) which is used to estimate the turbu¬ 
lent velocity of Perseus by Davis et al. (2008). 
Finally, using the cloud mass and turbnlent ve¬ 
locity we can obtain the turbulent momentum of 
~ 16 OOM 0 kms“^ in our survey region. 

Since the momentum of an outflow calculated 
using only a single emission line at 2.12 /rm is un¬ 
likely to be accurate, we also adopt canonical val¬ 
ues suggested by Davis et al. (2008) and loannidis 
& Froebrich (2012b) to estimate the momentum 
injected by these outflows. For typical magnitudes 
of outflow momentum in the range ~ 10 “^ and 
~ 2 X 10“® Mq kms“^ yr”^, for Class 0 and Class 
I sources, respectively (Bontemps et al. 1996), we 
estimate that a low-mass protostar will only in¬ 
ject ~1.0MQkms“^ (Davis et al. 2008; loanni¬ 
dis & Froebrich 2012b), over its lifetime that is 
usually on the order of 10^—10^ yrs. Our cal¬ 
culations therefore suggest, the sample of out¬ 
flows detected in Aquila can only possibly supply 
~ 4 OM 0 kms“^, which is a factor of 40 smaller 
than the observationally estimated cloud turbu¬ 
lent momentum. Evidently, the molecular hydro¬ 
gen outflows appear to be inefficient sources of tur¬ 
bulent momentum for this cloud. 

However, a caveat must be added here; the 
canonical values for outflow momentum suggested 
by Bontemps et al. (1996) were originally calcu¬ 
lated using measurements of “classical” slow CO 
molecular outflows. These suggested values of mo¬ 
mentum therefore only represent the low-velocity 
components of mass outflows. However, a large 
fraction of the jet momentum may be carried by 
a high-velocity collimated jet (Walawender et al. 
2005). It is therefore likely that we are underesti¬ 
mating the momentum injected by the MHOs in 
the Aquila molecular cloud. Nevertheless, it has 
been argued that the transport of the fast jet mo¬ 
mentum to the ambient cloud may be inefficient 


because the jets may be relatively unhindered by 
their surroundings, especially considering the exis¬ 
tence of parsec-scale jets and the high proper mo¬ 
tions of the distant jet knots (Davis et al. 2008). 

One possible way to balance the observed tur¬ 
bulent momentum of the molecular cloud is via 
many generations of outflows, as suggested by 
Davis et al. (2008) in Perseus. The lifetime of the 
giant molecular clouds is ~10^yr (Larson 1981), 
which is one or two orders of magnitude longer 
than the lifetime of Class O/I phase. Thus a few 
tens of generations of mass outflows may con¬ 
tribute to the observed turbulent momentum of 
the molecular clouds. 

5. Summary 

We have performed an unbiased near-infrared 
survey toward the Aquila molecular cloud in J, 
H, Kg, and H 2 bands using WIRCam equipped 
on CFHT. Our survey covers a sky region of ~ 
1° X1° which embraces the Serpens South cluster 
and W40 H II region. Combining the archival data 
from Spitzer and Herschel space telescopes, we 
investigate the properties of molecular hydrogen 
outflows in our survey region. The main results 
are listed as follows: 

1. We have identified 45 MHOs that consist of 
108 MHO features. Of 45 MHOs, 11 are 
previously known objects and 34 are newly 
discovered. Using the Spitzer archival data, 
we also identify 802 YSO candidates in our 
survey region. Based on the morphologies 
of MHOs and locations of MHOs and YSO 
candidates, we associate 43 MHOs with 40 
YSO candidates and obtain 40 molecular hy¬ 
drogen outflows. 

2. We obtain the visual extinction map of our 
survey region using the Herschel archival 
data through htting the SED from 160- 
500 /im pixel by pixel. Based on this extinc¬ 
tion map, we have tried to correct the fluxes 
of YSO candidates for extinction and clas¬ 
sify 802 YSO candidates into three classes 
with the de-reddened spectral indices: 100 
Class I sources, 551 Class H sources, and 
151 Class HI/MS sources. We find that 27 
(~ 68 %) of 40 outflow sources are protostars 
while ~27% of Class I sources drive H 2 out- 
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flows. We also find that 18 outflow sources 
are associated with infrared nebulae. 

3. In our outflow sample, the distribution 
of jet lengths shows a steep decrease in 
the number of outflows with increasing 
lengths, roughly following the relation of 

N OC I0{-0-S5±0.1)xlength{pc) ^ g^^pg 

of —0.85 is similar to the value of —0.75 
that is obtained with outflows from the 
low/intermediate mass stars in the Galactic 
plane toward Serpens and Aquila by loan- 
nidis & Froebrich (2012b). The distribution 
of jet position angles shows a nearly uniform 
distribution, which indicates that the molec¬ 
ular hydrogen outflows in the Aquila molec¬ 
ular cloud may orientate randomly. We also 
find that there is no obvious correlation be¬ 
tween jet lengths or jet opening angles and 
spectral indices of driving sources, which 
agrees with the results in Perseus (Davis 
et al. 2008), Orion A (Davis et al. 2009), 
and Ophiuchus (Zhang et al. 2013). 

4. We perform the areal photometry for our 
MHO features and estimate the H 2 1 — 
0 ^(l) line luminosities for the 40 outflows 
in the Aquila molecular cloud. The H 2 1 — 
0 >S'( 1 ) luminosities of our outflows range 
from ^3x10“® to ^ 2 . 6 x 10 “^ Lq with the 
median value of 4.7x 10-5 Lq, which is one 
or two order of magnitude lower than that in 
Orion A (Stanke et al. 2002) and that in the 
Galactic plane toward Serpens and Aquila 
(loannidis & Froebrich 2012a,b). We also 
find no obvious correlation between outflow 
H 2 1-0 ^(l) luminosities and spectral in¬ 
dices of driving sources. 

5. We estimate that the turbulent momen¬ 
tum in our survey region is about 1600 
M 0 kms-^. However, the momentum in¬ 
jection from molecular hydrogen outflows 
in our survey region is only ~40 Mq km s” ^, 
which is a factor of 40 smaller than the cloud 
turbulent momentum. Thus the momentum 
injection from molecular hydrogen outflows 
is not enough to match the observed turbu¬ 
lent momentum in our survey region of the 
Aquila molecular cloud. 
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Table 1 

Molecular hydrogen outelows identieied in Aquila. 





18 









o 

o 

o 


CM 


U 

0 ) 

o 


00 

• * 



- 05 ° 00’ 


18^ 34^ 32 30 28 

RA (J2000) 


Fig. 1.— Coverage of our data. The underlying image is the extinction map from Dobashi (2011). The blue 
contours represent the visual extinction of 5 and 10 magnitudes. Ten fields observed with CFHT/WlRCam 
are marked with red boxes while the black dashed box shows the coverage of our retrieved Spitzer archive 
data. 
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Fig. 2.— Left: Three-color image of the Aquila molecular cloud that is constructed with IRAC 3.6 /rm (blue), 
4.5/rm (green), and 8.0/rm (red) images. The YSOs identified in this paper are also marked with different 
symbols, red pentagrams indicating Class I sources, green triangles indicating Class II sources and cyan 
circles indicating transition disks. Right: Extinction map derived from Herschel PACS/SPIRE archive data 
of the Aquila molecular cloud with the resolution of 36.9", corresponding to the HPBW resolution of SPIRE 
at 500/im. The blue contours represent the visual extinction {Ay) of 4, 10, and 30 magnitudes. The location 
of the W40 H II region is marked with the white circle whose scale is in proportion to the size of W40 (Mallick 
et al. 2013). 
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Fig. 3.— Spatial distribution of molecular hydrogen outflows detected in this paper. The underlying is the 
Herschel extinction map obtained in Section 2.4. The 108 MHO features are marked with green circles and 
40 outflow sources are labeled with red filled pentagrams. The blue solid lines connect the outflow sources 
and their associated MHOs, hence indicating the identified molecular hydrogen outflows. The location of 
W40 H II region is marked with a black circle whose scale is in proportion to the size of W40 (Mallick et al. 
2013). 
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Fig. 4.— Left: Color-color diagrams showing all classified YSO candidates and field sources. The Class I 
candidates are marked with red triangles while the Class II candidates are labeled with green squares as well 
as the transition disk candidates are marked with blue pluses. The grey dots represent the Class Ill/field 
sources. Right: The H — Kg vs. J — H color-color diagram for the YSO candidates in the Aquila molecular 
cloud. The solid line show the intrincsic colors for the main-sequence stars (black) and giants (red) (Bessell 
& Brett 1988) while the dash-dotted line is the locus of T Tauri stars from Meyer et al. (1997). The dashed 
lines show the reddening direction and the arrow shows the reddening vector. The extinction law we adopted 
is from Wang & Jiang (2014). Note that the dashed lines separate the diagram into three regions marked 
with numbers 1, 2, and 3 in the figure. In each region, we use different methods to estimate the extinction 
of the YSO candidates (see the text for details). 
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Fig. 5.— Histograms of the apparent spectral indices (left panel) and de-reddened spectral indices (right 
panel) of all YSOs (hatched columns) and H 2 flow sources (filled black columns) identified in this paper. 
The dashed lines show the criteria of YSO classification scheme suggested by Lada (1987). 



Fig. 6.— Infrared nebulae associated with outflow sources. The backgrounds are iFs—band images. North 
is up and east to the left. The outflow sources are located in the center of panels. 
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Fig. 7.— A disk-jet system candidate that consists of MHO 3251bl and its driving source ID^17. The left 
panel shows the H 2 narrow-band image while the right panel shows the Kg broadband image. The white 
plus in the right panel marks the flux-weighted mean IRAC and MIPS 24/xm position for ID^17. The 
blue cross in the right panel marks the mean IRAC position for ID^17 while the red cross labels the MIPS 
24/im position for ID^I7. The solid line in the left panel shows the orientation of a collimated jet. 
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Fig. 8. — The region of MHO 3260. Top-left: three-color image constructed with CFHT J (blue), H 2 (green), 
and Ks (red) images. Top-right: three-color image constructed with IRAC 3.6/im (blue), 4.5 /rm (green), 
and 8.0 /rm (red) images. Bottom-left: three-color image constructed with MIPS 24 /rm (blue), PACS 
70 pm (green), and 160 /rm (red) images. Bottom-right: three-color image that is constructed with SPIRE 
250/im (blue), 350 /rm (green), and 500/im (red) images. The plus marks the position of Class 0 source, 
ID#2. 
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Fig. 9.— Distributions of jet lengths in bins of 0.2 pc {top left), jet opening angles in bins of 20° {top right), 
jet position angles in bins of 30° {bottom left), and logarithm of jet H 2 1 — 0 S'(l) luminosity in bins of 0.5 
{bottom right). 
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Fig. 10.— Distribution of jet lengths. The black filled circles show the distribution of jet lengths in bins 
of 0.15 pc while the solid line shows the linear fitting for this distribution. The gray dotted lines show the 
linear fittings for the distributions of jet lengths in different bins (0.1-0.3 pc in step of 0.01 pc) that have 
been scaled to the fitting in bins of 0.15 pc. 
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Fig. 11.— Jet length {left), jet opening angle {middle), and logarithm of jet H 2 1 — 0 5(1) luminosity 
{right) plot against source apparent spectral index {top) and source de-reddened spectral index {bottom). 
Note that outflow source ID#2 is invisible at 2-24 p,m and ID^19 is saturated in all IRAC bands and MIPS 
24^m band. To plot these sources on the figure, we assume the value of 3 for their spectral indices because 
they are associated with the protostars identified by Bontemps et al. (2010). 
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A. Description of the H2 outflows in Aquila 

Finally we have identified 45 MHOs which consist of 108 MHO features in the Aquila molecular cloud. 
Table A1 shows the information of these MHO features, including their positions and areas. Note that we 
use a polygon to define the morphology of each MHO feature. The position and area of each MHO feature 
in table A1 are calculated based on the polygon. We also present a continuum-subtracted H 2 image and 
give a brief description for each MHO. Figs. A1-A42 show these continuum-subtracted H 2 images. The blue 
and pink polygons mark the MHO features in the figures. Note that for some MHO features there are holes 
that are labeled with red circles in the polygons. The MHO features which belong to the same H 2 outflow 
and their possible driving source are connected with black dashed lines in the figures. 

• MHO 2213: Fig. A1 shows the region of MHO 2213 which consists of three features, MHO 2213a-c. 
Connelley et al. (2007) firstly detected this object, however, they only detected MHO 2213a due to 
the small sky coverage (~l'xl') of the observation. Connelley et al. (2007) also identihed an infrared 
nebula which is associated with IRAS 18264-0143 and visible in our Ag—band image (see Fig. 6). We 
have identihed a YSO, ID#3, which is also associated with IRAS 18264-0143. It seems that MHO 2213 
and ID#3 constitute a bipolar outflow system. Therefore, we suggest ID#3 as the possible driving 
source of MHO 2213. 

• MHO 2214: Connelley et al. (2007) flrstly detected this object, but only the feature that is labeled as 
MHO 2214b in Fig. A2. Teixeira et al. (2012) increased the number of MHO features up to two, MHO 
2214b and MHO 2214r, corresponding to MHO 2214b and MHO 2214r5 in Fig. A2, respectively. Using 
our deep H 2 images, we have identihed six features in MHO 2214, MHO 2214rl-r5 and MHO 2214b. 
Teixeira et al. (2012) suggested a protostar, “PO” (Bontemps et al. 2010), as the possible driving source 
of MHO 2214 because “PO” and MHO 2214 are nearly located on a line. We identify a YSO, ID#9, 
which is associated with “PO”. Thus we accept the conclusion of Teixeira et al. (2012) and suggest 
ID^9 as the possible driving source of MHO 2214. 

• MHO 3247: Fig. A3 shows the region of MHO 3247 that consists of three features, MHO 3247a- 
c. Teixeira et al. (2012) hrstly identihed MHO 3247, but they only detected the northeast lobe, 
corresponding to MHO 3247c in Fig. A3. Using our deep H 2 images, we detect the southwest lobe 
of MHO 3247 that is labeled as MHO 3247a and b. Teixeira et al. (2012) suggested the protostar, 
“P2”(Bontemps et al. 2010), as the possible driving source of MHO 3247. We also identify a YSO, 
ID^21, which is associated with “P2”. Thus we accept the conclusion of Teixeira et al. (2012) and 
also suggest ID#21 as the possible driving source of MHO 3247. 

• MHO 3248 and MHO 3251b2: Fig. A4 shows the region of MHO 3248 and MHO 3251b2. Teixeira 
et al. (2012) hrstly detected part of features of MHO 3248 which mainly corresponds to MHO 3248a 
in Fig. A4. They suggested the Class 0 source, “P3”, as the driving source of MHO 3248. We detect 
another lobe of MHO 3248, labeled as MHO 3248b-c. We also identify a YSO, ID^22, which is 
associated with “P3”. Thus we accept the suggestion of Teixeira et al. (2012) and suggest ID#22 as 
the possible driving source of MHO 3248. MHO 3251b2 is hrstly detected by Teixeira et al. (2012), 
too. Teixeira et al. (2012) suggested that MHO 3251b2 may result from the interaction of two different 
outhows, among which one consists of “P3” and MHO 3248 and the other consists of MHO 3251 (see 
the description about MHO 3251). MHO 3251b2 shows more detailed structures in our deep H 2 image. 
Based on the morphology of MHO 3251b2, we propose that MHO 3251b2 belongs to the outhow that 
is driven by ID=f^22. Thus we suggest that MHO 3248, MHO 3251b2, and ID^22 constitute a bipolar 
outhow system. 

• MHO 3249: Teixeira et al. (2012) hrstly detected this object and suggested the Class I source, 
“P4”(Bontemps et al. 2010), as the possible driving source. We also identify a YSO, ID#24 that 
is associated with “P4”. Thus here we suggest ID#24 as the driving source of MHO 3249, too. Note 
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that we just mark the heads of two lobes of the outflow that is driven by ID^24 in Fig. A5 with the 
blue polygons. In fact, Teixeira et al. (2012) suggested that there may be interaction between MHO 
3248 and the blueshifted lobe of MHO 3249. They also suggested that part of MHO 3248a may belong 
to the outflow that consists of MHO 3249 and ID#24. MHO 3248a shows more detailed structures in 
our deep H 2 image and obviously consists of two components in Fig. A4 and Fig. A5: a north-south 
component may belong to the outflow driven by ID#22; a NE-SW component should belong to the 
outflow driven by ID^24. However, these two components are overlapped and hard to distinguish from 
each other. Thus we do not mark the sub-structures of MHO 3248a in Fig. A4. Note especially that 
we assume that the whole MHO 3248a belongs to the outflow driven by ID^22 when we calculate the 
outflow parameters such as jet length, jet opening angle, and jet position angle. 

• MHO 3250: Teixeira et al. (2012) firstly detected this object, but only the feature that is marked as 
MHO 3250c in Fig. A6. We identify more features, including MHO 3250a, b and d. Teixeira et al. 
(2012) suggested a young star, “Yl” (Bontemps et al. 2010), as the driving source of MHO 3250. We 
have identified a YSO, ID^18, which is associated with “Yl”. Thus here we accept the conclusion of 
Teixeira et al. (2012) and suggest that MHO 3250 is driven by ID:j(^:18. 

• MHO 3251: Fig. A7 shows the region of MHO 3251. Teixeira et al. (2012) firstly detected this object 
and suggested the Class I source, “Y2” (Bontemps et al. 2010) that is associated with our identified 
YSO ID#17, as the possible driving source. Teixeira et al. (2012) suggested that the redshifted lobe of 
the outflow driven by ID^17 could extend to MHO 3250c (see Fig. A6) and there may be interaction 
between the outflow driven by ID^17 and the outflow driven by ID^18 (see the description about 
MHO 3250). MHO 3250c shows more detailed structures in our deep H 2 image (see Fig. A6) and 
we argue that there is no obvious indication of interaction based on the morphology of MHO 3250c. 
Moreover, we detect a new feature, MHO 3251r, being located between ID^17 and MHO 3250c. Thus 
we suggest MHO 3251r as the outskirt of the redshifted lobe of the outflow driven by ID^17. On the 
other hand, Teixeira et al. (2012) suggested that the blueshifted lobe of the outflow driven by ID#17 
could extend to MHO 3251b2. However, we argue that MHO 3251b2 is likely to belong to the outflow 
driven by ID#22 (see the description about MHO 3248). We detect a new feature, MHO 3251bl that 
shows a highly collimated jet-like structure between ID^17 and MHO 3251b2 in Fig. A7. Therefore, 
we suggest MHO 3251bl as the outskirt of the blueshifted lobe of the outflow driven by ID^17. 

• MHO 3252: Teixeira et al. (2012) firstly detected this object, but only the features that correspond to 
MHO 3252bl-b2, and r4 in Fig. A8. We identify three new features in the west of MHO 3252r4, MHO 
3252rl-r3. Teixeira et al. (2012) suggested the protostar, “Y3” (Bontemps et al. 2010), as the possible 
driving source of MHO 3252. We detected an infrared source, ID^19, as the counterpart of protostar 
“Y3”. ID^19 is visible at 2-24/xm, but it is saturated in all IRAC and MIPS 24 /im bands. Based on 
the good alignment of MHO 3252 and ID#19, we suggest ID#19 as the driving source of MHO 3252. 

• MHO 3253: Fig. A9 shows the region of MHO 3253. Teixeira et al. (2012) firstly detected this object, 
but only the feature corresponding to MHO 3253a in Fig. A9. We detect a new feature, MHO 3253b, 
in the east of MHO 3253a. Teixeira et al. (2012) suggested the Class I source “Y4” (Bontemps et al. 
2010), which is associated with our identified YSO ID^28, as the driving source of MHO 3253. Here 
we accept the suggestion of Teixeira et al. (2012) and propose ID#28 as the driving source of MHO 
3253. 

• MHO 3254 and MHO 3277-3278: Fig. A9 also shows the region of MHO 3254 and MHO 3277-3278. 
MHO 3254 was firstly identified by Teixeira et al. (2012), but Teixeira et al. (2012) only detected two 
features that correspond to MHO 3254c in Fig. A9 and MHO 3254a in Fig. AlO. We detect a new 
feature, MHO 3254b, lying between MHO 3254a and c. Teixeira et al. (2012) suggested a YSO that 
is located to the east of ID#28, “Y5”, as the possible driving source of MHO 3254a and c. We argue 
that MHO 3254c shows bow-like structure in our deep H 2 image and seems to be able to constitute. 
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together with MHO 3253 and ID^28, a bipolar outflow system. We also detect two new MHOs in the 
southwest of ID#28, MHO 3277 and MHO 3278. It seems that MHO 3253, ID#28, MHO 3254b-c, 
and MHO 3277-3278 align well with each other. Thus here we propose another possibility that MHO 
3253, MHO 3254b-c, and MHO 3277-3278 are driven by ID:j^28. Moreover, based on the location of 
MHO 3254a, we propose that MHO 3254a belongs to the outflow driven by ID#23 (see the description 
about MHO 3255). 

• MHO 3255: Teixeira et al. (2012) firstly detected this object, but only the features corresponding to 

MHO 3255a and b in Fig. AlO. We identify a new feature, MHO 3255c, which is located ~4' to the 
northeast of MHO 3255b. Teixeira et al. (2012) suggested a protostar “P5” that is associated with our 
identified YSO, ID#23, as the possible driving source of MHO 3255. Based on the good alignment 
among MHO 3255, MHO 3254a, and we suggest that ID:j(^:23 drives a bipolar outflow that 

consists of MHO 3255 and MHO 3254a. 

• MHO 3259: Fig. All shows the region of MHO 3259 that exhibits a bow-like structure. ID^l is the 
nearest YSO, lying in the southwest of MHO 3259. We propose ID^l as the possible driving source 
of MHO 3259. 

• MHO 3260: Fig. A12 shows the region of MHO 3260 which consists of four features. MHO 3260c and 
d are two knots. MHO 3260a and b show approximate point symmetric structures, which indicates 
that the driving source should be located between them. Bontemps et al. (2010) detected a protostar 
which is labeled as ID^2 in Fig. A12. Note that ID^2 is invisible at 3.6-24/im, but visible in Herschel 
PACS and SPIRE images, which indicates that ID#2 may stay at a very early evolutionary stage. We 
suggest ID#2 as the driving source of MHO 3260. 

• MHO 3261: Fig. A13 shows the region of MHO 3261 which consists of two features. The YSO ID^4 
is located at the northeast of MHO 3261b and associated with an infrared nebula (see Fig. 6). The 
infrared nebula exhibits the NE-SW elongated structure. We suggest ID#4 as the driving source of 
MHO 3261. 

• MHO 3262: This object is an elongated knot in Fig. A14. ID#5 is the nearest YSO to it and we 
suggest ID#5 as the possible driving source of MHO 3262. 

• MHO 3263: Fig. A15 shows the region of MHO 3263 which consists of three features. The YSO ID^6 
and MHO 3263 are roughly located on a straight line. Thus we suggest ID^6 as the possible driving 
source of MHO 3263. 

• MHO 3264: Fig. A16 shows the region of MHO 3264 which consists of five features. The YSO ID^7 is 
roughly aligned with the features of MHO 3264. ID#7 is also associated with an infrared nebula (see 
Fig. 6) which seems to trace the cavity wall of the west lobe of outflow. Thus we suggest ID#7 as the 
possible driving source of MHO 3264. 

• MHO 3265: Fig. A17 shows the region of MHO 3265 which consists of two features. The YSO ID^8 
is located between MHO 3265a and b. Note that we suggest MHO 3265a belongs to another outflow 
rather than the outflow that consists of MHO 2214 based on its morphology although MHO 3265a is 
located between MHO 2214r3 and MHO 2214r4 (see the description about MHO 2214). Considering 
the good alignment of ID:j(^:8 and MHO 3265, we propose ID^8 as the possible driving source of MHO 
3265. 

• MHO 3266: MHO 3266 shows a complex structure in Fig. A18. The YSO ID:j(^:10 is located in the 
southwest of MHO 3266 and seems to be located on the extension line of the major axis of MHO 3266. 
Thus we suggest ID#10 as the possible driving source of MHO 3266. 
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• MHO 3267: Fig. A19 shows the region of MHO 3267 which consists of three features. The YSO 

is associated with IRAS 18270-0153 (Connelley et al. 2007; Connelley & Greene 2010) and an infrared 
nebula (Connelley et al. 2007). This infrared nebula shows a NE-SW bipolar structure in our Kg image 
(see Fig. 6). It seems that ID^ll, the infrared nebula, and MHO 3267 constitute a bipolar outflow 
system. 

• MHO 3268: This object consists of four features, MHO 3268a-d, in an arrangement of a fan shape with 

the center of in Fig. A20. is also associated with a protostar identified by Bontemps 

et al. (2010) using the Herschel data. Here we suggest ID#12 as the possible driving source of MHO 
3268. 

• MHO 3269: Fig. A21 shows the region of MHO 3269. The YSO ID:;(^13 is located to the east of MHO 
3269 and associated with an infrared nebula that seems to trace the cavity walls of an east-west outflow 
(see Fig. 6). MHO 3269 and ID^13 are roughly on a straight line. Thus we suggest ID^13 as the 
possible driving source of MHO 3269. 

• MHO 3270: Fig. A22 shows the region of MHO 3270 that consists of two features. MHO 3270a seems 

to consist of a knot and a jet-like structure driven by the YSO while MHO 3270b shows a 

fan-shaped structure. ID^14 and MHO 3270 align well with each other. Thus we suggest ID^14 as 
the possible driving source of MHO 3270. 

• MHO 3271: Fig. A23 shows the region of MHO 3271. There are several YSOs in the ambient of MHO 
3271. Among them, ID^15 is associated with an infrared nebula (see Fig. 6). It seems that ID:;(^:15, the 
infrared nebula, and MHO 3271 constitute a NE-SW bipolar outflow system and the infrared nebula 
represents the northeast lobe while MHO 3271 traces the southwest lobe. 

• MHO 3272: MHO 3272 shows a complex structure in Fig. A24. ID^16 is the nearest YSO. It seems 
that ID^16 and MHO 3272 constitute a bipolar outflow system. 

• MHO 3273: Fig. A25 shows the region of MHO 3273 which consists of eight features. The YSO 1T>^2{) 
is located in the vicinity of MHO 3273d. ID41^2Q and MHO 3273 are located roughly on a straight line. 
Based on the morphologies and locations, we suggest ID^20 as the possible driving source of MHO 
3273. 

• MHO 3274: Fig. A26 shows the region of MHO 3274. ID#25 is the nearest YSO and associated with 

an infrared nebula which is shown in Fig. 6. It seems that lT)4l^2b^ the infrared nebula, and MHO 3274 
constitute a bipolar outflow system. Thus we suggest as the driving source of MHO 3274. 

• MHO 3275: This object is an elongated knot in Fig. A27. lT>^2b (see Fig. A26) is the nearest YSO, 
but ID^25 drives an east-west outflow (see the description about MHO 3274), which disagrees with 
the location of MHO 3275. Therefore, we suggest the secondary nearest YSO, ID#26, as the possible 
driving source of MHO 3275. 

• MHO 3276: This object exhibits the bow-like structure in Fig. A28, which indicates that its possible 
driving source may be located to the southwest of it. ID^27 is located in the southwest of MHO 3276 
and also associated with a protostar identified by Bontemps et al. (2010) using the Herschel data. Thus 
we suggest ID#27 as the possible driving source of MHO 3276. 

• MHO 3279: Fig. A29 shows the region of MHO 3279 that consists of two features. The YSO 1T>^2Q 
is located between the two features of MHO 3279. It seems that MHO 3279 and ID4i^2^ constitute a 
bipolar outflow system. Note that ID^29 is also associated with a SE-NW infrared nebula (see Fig. 6). 
Thus we suggest ID#29 as the driving source of MHO 3279. 
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• MHO 3280: This object shows bow shock morphology in Fig. A30, which indicates that its possible 
driving source may be located to its northwest. The YSO ID^30 is located in the northwest of MHO 
3280 and associated with an infrared nebula that exhibits SE-NW bipolar structure in Fig. 6. We 
suggest ID#30 as the driving source of MHO 3280. 

• MHO 3281: Fig. A31 shows the region of MHO 3281 that consists of three features. The YSO ID:j^31 
is located between MHO 3281a and b. Note that ID^31 is also associated with an infrared nebula (see 
Fig. 6). It seems that MHO 3281 and ID:j(131 constitute a bipolar outflow system. 

• MHO 3282: This object seems to be a bright knot with a faint ring-like structure in Fig. A32. 
is the nearest YSO and we suggest it as the possible driving source of MHO 3282. 

• MHO 3283: Fig. A33 shows the region of MHO 3283 that consists of two features. MHO 3283a exhibits 
the multiple bow shock morphology, which indicates that its possible driving source should be located 
to its east. The YSO ID^33 in the southeast of MHO 3283 is located roughly on the line connecting 
MHO 3283a and MHO 3283b. Thus we suggest ID^33 as the possible driving source of MHO 3283. 

• MHO 3284 and MHO 3285: MHO 3284 shows complex morphology in Fig. A34 while MHO 3285 
consists of some knots in Fig. A35. No YSOs have been identified in the ambient of MHO 3284 and 
MHO 3285. Therefore, we can not infer any driving source candidates for MHO 3284 and MHO 3285. 

• MHO 3286: This object is a chain of two bright knots in Fig. A36, which indicates that its possible 
driving source may be located to its southeast or northwest. There is a YSO, ID#34, to the southeast 
of MHO 3286. ID^34 and the knots of MHO 3286 seems to align with each other. Thus we suggest 
ID#34 as the possible driving source of MHO 3286. 

• MHO 3287: Figure A37 shows the region of MHO 3287 that consists of three features. MHO 3287a 
is a collimated jet in the vicinity of the YSO ID^35 that is associated with a protostar identified by 
Bontemps et al. (2010) using the Herschel data. There is good alignment between ID^35 and MHO 
3287. Thus we suggest ID#35 as the driving source of MHO 3287. Note that the circular bright blot 
in the corresponding position of ID#35 is due to the instrument effect, “persistence”, which occurs 
due to charges being trapped in high flux or saturated regions in one exposure that are then slowly 
released in subsequent exposures. Since exposures are dithered, the released charges appear as new 
“artificial sources” in the subsequent images. 

• MHO 3288: This object is a patch in Fig. A38. The nearest YSO to MHO 3288 is ID^36, which is 
also associated with an infrared nebula (see Fig. 6). We suggest ID#36 as the possible driving source 
of MHO 3288. 

• MHO 3289: This object shows a bow-like structure in Fig. A39, which indicates that the driving source 
of MHO 3289 may be located to the southeast of it. ID#37 is located to the southeast of MHO 3289 
and associated with a SE-NW bipolar infrared nebula (see Fig. 6). We suggest ID#37 as the possible 
driving source of MHO 3289. 

• MHO 3290: This object consists of three features that align well with each other in Fig. A40. The 
YSO ID^38 is located between MHO 3290b and MHO 3290c and associated with an infrared nebula 
that extends towards southeast from ID#38 in Fig. 6. It seems that ID^38 and MHO 3290 constitute 
a bipolar outflow system. We suggest ID#38 as the driving source of MHO 3290. 

• MHO 3291: This object is a knot in Fig. A41. ID^39 is the nearest YSO to MHO 3291. We suggest 
ID^39 as the possible driving source of MHO 3291. 

• MHO 3292: Figure A42 shows the region of MHO 3292. The YSO is located roughly on the line 

connecting the features of MHO 3292. ID^40 is also associated with an infrared nebula that extends 
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towards the northeast from ID^40 in Fig. 6. Thus we suggest as the possible driving source of 

MHO 3292. 
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Table A1 

MHO FEATURES IN AqUILA. 


MHO 

RA"* 

(J2000) 

Dec^ 

(J2000) 

Area^ 

(arcsec^) 

Radius*^ 

(arcsec) 

H 2 1 - 0 S(l) flux 
(10-i» W m"^) 

3259 

18 

28 

41.4 

-01 

32 

08 

24.1 

2.8 

19.31±1.85 

3261a 

18 

28 

52.8 

-01 

46 

48 

9.5 

1.7 

2.11±0.68 

3260a 

18 

29 

03.0 

-01 

39 

10 

31.2 

3.2 

10.57±0.89 

3260b 

18 

29 

04.4 

-01 

39 

01 

33.6 

3.3 

27.72T1.98 

2213a 

18 

29 

04.9 

-01 

41 

52 

61.9 

4.4 

33.70±2.61 

3260c 

18 

29 

05.4 

-01 

38 

52 

5.8 

1.4 

0.69±0.13 

3260d 

18 

29 

06.2 

-01 

38 

47 

5.4 

1.3 

l.OOitO.16 

2213b 

18 

29 

07.9 

-01 

42 

16 

26.5 

2.9 

6.49±0.64 

2213c 

18 

29 

08.8 

-01 

42 

34 

57.8 

4.3 

29.99±2.35 

3261b 

18 

29 

08.9 

-01 

46 

27 

34.5 

3.3 

4.83±1.48 

3262 

18 

29 

14.1 

-01 

57 

59 

13.3 

2.1 

1.87±0.62 

3252rl 

18 

29 

15.8 

-02 

13 

21 

44.3 

3.8 

10.08T1.75 

3264a 

18 

29 

17.8 

-01 

39 

17 

50.0 

4.0 

16.65±3.96 

3268a 

18 

29 

18.9 

-02 

19 

37 

370.1 

10.9 

11.30±2.34 

3264b 

18 

29 

20.6 

-01 

38 

59 

102.7 

5.7 

10.63±2.67 

3263a 

18 

29 

21.2 

-01 

37 

33 

9.1 

1.7 

1.23±0.37 

3263b 

18 

29 

22.5 

-01 

37 

52 

32.5 

3.2 

3.71±1.01 

3268b 

18 

29 

24.9 

-02 

19 

36 

571.3 

13.5 

75.38±11.76 

2214rl 

18 

29 

25.6 

-01 

49 

55 

40.8 

3.6 

2.18±0.41 

3264c 

18 

29 

25.7 

-01 

38 

59 

28.8 

3.0 

2.75±0.80 

3268c 

18 

29 

26.5 

-02 

18 

28 

254.8 

9.0 

219.70±32.37 

3263c 

18 

29 

26.9 

-01 

39 

49 

17.6 

2.4 

3.35±0.90 

2214r2 

18 

29 

27.8 

-01 

50 

03 

61.5 

4.4 

5.26±0.78 

3264d 

18 

29 

29.3 

-01 

38 

41 

39.3 

3.5 

6.80±1.72 

3252r2 

18 

29 

30.1 

-02 

12 

20 

121.1 

6.2 

59.08±9.05 

3268d 

18 

29 

30.5 

-02 

16 

27 

59.0 

4.3 

8.69±1.64 

3273a 

18 

29 

32.5 

-01 

32 

35 

251.1 

8.9 

64.67±14.59 

3252r3 

18 

29 

33.0 

-02 

11 

52 

103.1 

5.7 

20.26±3.34 

2214r3 

18 

29 

33.3 

-01 

50 

33 

81.2 

5.1 

12.67T1.59 

3269a 

18 

29 

33.6 

-01 

33 

12 

34.5 

3.3 

5.91±1.52 

3265a 

18 

29 

34.0 

-01 

50 

28 

25.5 

2.8 

8.35±1.08 

2214r4 

18 

29 

34.8 

-01 

50 

31 

43.8 

3.7 

15.19±1.83 

3267a 

18 

29 

36.4 

-01 

52 

00 

5.0 

1.3 

0.67±0.14 

2214r5 

18 

29 

36.4 

-01 

50 

43 

129.9 

6.4 

38.15±4.27 

3267b 

18 

29 

36.9 

-01 

51 

38 

8.4 

1.6 

2.48±0.38 

3269b 

18 

29 

38.2 

-01 

33 

02 

36.7 

3.4 

16.19±3.81 

2214b 

18 

29 

39.7 

-01 

51 

27 

186.1 

7.7 

108.30±11.38 

3265b 

18 

29 

39.7 

-01 

52 

46 

312.6 

10.0 

10.11±1.53 

3267c 

18 

29 

40.1 

-01 

51 

08 

46.0 

3.8 

5.75±0.82 

3269c 

18 

29 

40.3 

-01 

32 

58 

19.4 

2.5 

4.75±1.23 

3250a 

18 

29 

43.2 

-02 

05 

29 

4.9 

1.2 

1.74±0.36 

3264e 

18 

29 

45.0 

-01 

37 

36 

18.0 

2.4 

9.69±2.35 

3270b 

18 

29 

47.0 

-01 

59 

34 

57.8 

4.3 

18.95±1.20 

3270a 

18 

29 

47.0 

-01 

55 

52 

8.9 

1.7 

2.73±0.26 

3250b 

18 

29 

50.8 

-02 

05 

08 

58.1 

4.3 

4.75±0.93 

3273b 

18 

29 

50.9 

-01 

34 

38 

294.3 

9.7 

43.96±10.16 

3266 

18 

29 

53.0 

-01 

29 

27 

436.9 

11.8 

204.70±55.70 

3273c 

18 

29 

55.3 

-01 

35 

38 

102.2 

5.7 

17.63±4.18 

3250c 

18 

29 

56.8 

-02 

04 

33 

333.1 

10.3 

95.79±14.45 

3271 

18 

29 

57.3 

-02 

02 

28 

44.4 

3.8 

3.30±0.66 

3278 

18 

29 

57.8 

-02 

12 

55 

32.0 

3.2 

7.62±1.35 

3252r4 

18 

29 

58.2 

-02 

10 

52 

243.0 

8.8 

70.21±10.72 

3272a 

18 

29 

59.2 

-02 

25 

53 

216.5 

8.3 

74.74±10.09 

3247a 

18 

30 

00.0 

-02 

03 

12 

23.0 

2.7 

9.79±1.66 

3251r 

18 

30 

00.1 

-02 

05 

26 

102.6 

5.7 

11.29±1.98 

3272b 

18 

30 

00.7 

-02 

24 

08 

137.8 

6.6 

44.77±6.16 

3249b 

18 

30 

00.8 

-02 

04 

59 

83.5 

5.2 

73.10±11.07 

3252bl 

18 

30 

01.8 

-02 

10 

24 

29.0 

3.0 

50.41±7.73 
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Table A1— Continued 


MHO 

RA"* 

(J2000) 

Dec^ 

(J2000) 

Area^ 

(arcsec^) 

Radius*^ 

(arcsec) 

H 2 1 - 0 S(l) flux 
(10-18 W m-^) 

3251bl 

18 

30 

02.0 

-02 

06 

28 

7.5 

1.5 

1.06±0.26 

3273d 

18 

30 

02.8 

-01 

36 

30 

30.7 

3.1 

8.16±2.03 

3247b 

18 

30 

03.0 

-02 

02 

59 

123.6 

6.3 

192.30±28.32 

3277 

18 

30 

03.1 

-02 

12 

18 

94.5 

5.5 

9.72±1.74 

3248a 

18 

30 

03.5 

-02 

03 

57 

1782.2 

23.8 

807.50±116.70 

3248b 

18 

30 

03.8 

-02 

01 

48 

114.4 

6.0 

20.75±3.36 

3255a 

18 

30 

03.9 

-02 

11 

11 

50.4 

4.0 

16.18±2.68 

3247c 

18 

30 

04.2 

-02 

02 

26 

225.2 

8.5 

155.50±23.05 

3248c 

18 

30 

04.2 

-02 

02 

05 

164.3 

7.2 

24.79±3.99 

3255b 

18 

30 

05.7 

-02 

10 

52 

308.7 

9.9 

308.80±45.18 

3251b2 

18 

30 

05.8 

-02 

08 

03 

1799.0 

23.9 

372.00±54.64 

3249r 

18 

30 

07.8 

-02 

01 

04 

64.8 

4.5 

12.56±2.14 

3254a 

18 

30 

08.3 

-02 

10 

49 

62.3 

4.5 

21.92±3.54 

3273e 

18 

30 

08.7 

-01 

37 

44 

18.8 

2.4 

2.44±0.71 

3273f 

18 

30 

11.1 

-01 

38 

09 

17.6 

2.4 

1.34±0.44 

3252b2 

18 

30 

14.1 

-02 

09 

22 

2816.7 

29.9 

316.00±46.85 

3250d 

18 

30 

16.1 

-02 

01 

52 

294.1 

9.7 

25.90±4.21 

3274 

18 

30 

16.4 

-02 

07 

27 

28.2 

3.0 

1.98±0.46 

3275 

18 

30 

16.7 

-02 

06 

28 

10.4 

1.8 

1.40±0.31 

3254b 

18 

30 

17.8 

-02 

10 

58 

128.1 

6.4 

20.26±3.37 

3255c 

18 

30 

21.8 

-02 

09 

56 

68.8 

4.7 

5.36±1.07 

3254c 

18 

30 

22.6 

-02 

10 

57 

108.8 

5.9 

33.86±5.33 

3253a 

18 

30 

27.2 

-02 

10 

35 

26.6 

2.9 

18.44±2.98 

3253b 

18 

30 

28.1 

-02 

10 

26 

33.3 

3.3 

2.45±0.54 

3279a 

18 

30 

28.3 

-01 

55 

59 

10.8 

1.9 

3.37±0.95 

3276 

18 

30 

28.9 

-01 

48 

19 

37.5 

3.5 

13.63±3.21 

3279b 

18 

30 

29.1 

-01 

56 

07 

22.8 

2.7 

2.11±0.65 

3273g 

18 

30 

43.9 

-01 

42 

05 

106.8 

5.8 

14.81±3.60 

3273h 

18 

30 

46.5 

-01 

43 

29 

28.2 

3.0 

7.61±1.90 

3283a 

18 

30 

46.9 

-02 

03 

23 

578.8 

13.6 

55.04±10.94 

3281a 

18 

30 

48.2 

-01 

56 

32 

16.9 

2.3 

6.35±0.84 

3281b 

18 

30 

49.0 

-01 

55 

45 

15.6 

2.2 

3.15±0.48 

3281c 

18 

30 

50.4 

-01 

55 

38 

41.4 

3.6 

5.61±0.79 

3283b 

18 

30 

51.8 

-02 

03 

52 

487.6 

12.5 

75.01±14.60 

3280 

18 

30 

55.1 

-01 

57 

37 

32.8 

3.2 

14.09±1.68 

3282 

18 

30 

59.4 

-02 

29 

19 

72.6 

4.8 

16.27±2.13 

3286 

18 

31 

21.7 

-02 

30 

12 

35.6 

3.4 

10.59±1.72 

3284 

18 

31 

22.0 

-02 

19 

58 

186.8 

7.7 

96.04±10.74 

3285 

18 

31 

24.6 

-02 

19 

55 

131.8 

6.5 

58.25±6.68 

3287a 

18 

31 

51.1 

-02 

29 

40 

80.6 

5.1 

6.71±1.11 

3287b 

18 

31 

58.9 

-02 

29 

22 

200.1 

8.0 

10.99±1.74 

3288 

18 

32 

05.1 

-02 

21 

18 

15.5 

2.2 

1.61±0.26 

3287c 

18 

32 

08.5 

-02 

29 

12 

211.0 

8.2 

17.84±2.55 

3289 

18 

32 

08.9 

-01 

55 

13 

29.9 

3.1 

10.18±0.88 

3290a 

18 

32 

11.5 

-02 

33 

16 

341.6 

10.4 

25.68±2.91 

3290b 

18 

32 

14.7 

-02 

34 

07 

595.7 

13.8 

38.09±4.16 

3291 

18 

32 

20.1 

-02 

34 

45 

10.0 

1.8 

3.22±0.44 

3290c 

18 

32 

20.2 

-02 

36 

05 

47.5 

3.9 

16.49±1.82 

3292a 

18 

32 

23.4 

-01 

52 

39 

35.8 

3.4 

19.24±1.87 

3292b 

18 

32 

28.7 

-01 

52 

19 

20.7 

2.6 

3.21±0.43 


®^Center of the polygon that is identifited based on the morphology and surface 
brightness distribution of each MHO feature. 

^Area of the polygon. 

‘^Equivalent circular radius of the polygon. 
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Fig. Al.— Region of MHO 2213. The background is the continuum-subtracted H 2 image and the MHO 
features are marked with blue and pink polygons. The black dashed line connects the MHO features that 
belongs to the same molecular hydrogen outflow and their possible driving source that is labeled with red 
filled pentagram. The number marked around the pentagram corresponds to the ID in column 1 of table 1. 
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Fig. A2.— Same as Fig. Al, but for the region of MHO 2214. Note that the square images titled with the 
names of MHO features show the zoomed views of MHO features which belong to the same H 2 outflow. 
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Fig. A3.— Same as Fig. Al, but for the region of MHO 3247. 
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Fig. A4.— Same as Fig. A2, but for the region of MHO 3248 and MHO 3251b2. Note that the red circles mark the holes inside polygons. 
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Fig. A5.— Same as Fig. A2, but for the region of MHO 3249. Note that the red circle marks the hole inside 
the polygon. 
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Fig. A6.— Same as Fig. A2, but for the region of MHO 3250. Note that the red circles mark the holes inside 
the polygons. 
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Fig. A7.— Same as Fig. Al, but for the region of MHO 3251. 
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Fig. A8.— Same as Fig. A2, but for the region of MHO 3252. Note that the red circles mark the holes inside 
the polygons. 


43 


















(oooir) (COOtr) MQ (ooozr) 3*0 


tM032V3<> 


»»«32Wo 





-oi* la so- V 

1«* ».»• J7-S* 

M (XKOO) 
»»4032&4b 


Iff «ff- 


-02* n- 00 ^ 



-ca* n- J<r 



!»" XT M.O* 2»ff* 
M (jmo) 
MM03277 


Iff *0-1 


nf* XT 24^ 

W (JKOO) 

MH0327B 


-oa* ir JO-I 



-02” 13- 30" 


18^ 30"^ 30* 


RA (J2000) 

Fig. A9.— Same as Fig. A2, but for the region of MHO 3253, MHO 3254b, and MHO 3277-3278. 
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Fig. AlO.— Same as Fig. A2, but for the region of MHO 3254a and MHO 3255. Note that the red circle 
marks the hole inside the polygon. 
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Fig. All.— Same as Fig. Al, but for the region of MHO 3259. 
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Fig. A12.— Same as Fig. A2, but for the region of MHO 3260. 
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Fig. A13.— Same as Fig. A2, but for the region of MHO 3261. 
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Fig. A14.— Same as Fig. Al, but for the region of MHO 3262. 
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Fig. A15.— Same as Fig. A2, but for the region of MHO 3263. 
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Fig. A17.— Same as Fig. A2, but for the region of MHO 3265. 
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Fig. A18.— Same as Fig. Al, but for the region of MHO 3266. Note that the red circles mark the holes 
inside the polygons. 
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Fig. A19.— Same as Fig. A2, but for the region of MHO 3267. 
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Fig. A21.— Same as Fig. A2, but for the region of MHO 3269. Note that the red circle marks the hole inside the polygon. 
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Fig. A22.— Same as Fig. A2, but for the region of MHO 3270. Note that the red circle marks the hole 
inside the polygon. 
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Fig. A24.— Same as Fig. Al, but for the region of MHO 3272. Note that the red circles mark the holes 
inside the polygons. 
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Fig. A25.— Same as Fig. A2, but for the region of MHO 3273. Note that the red circles mark the holes 
inside the polygons. 
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Fig. A26.— Same as Fig. Al, but for the region of MHO 3274. 
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Fig. A27.— Same as Fig. Al, but for the region of MHO 3275. 
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Fig. A28.— Same as Fig. A2, but for the region of MHO 3276. 
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Fig. A29.— Same as Fig. Al, but for the region of MHO 3279. 
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Fig. A30.— Same as Fig. Al, but for the region of MHO 3280. 
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Fig. A31.— Same as Fig. A2, but for the region of MHO 3281. 
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Fig. A33.— Same as Fig. A2, but for the region of MHO 3283. Note that the red circles mark the holes 
inside the polygons. 
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Fig. A34.— Same as Fig. Al, but for the region of MHO 3284. Note that we did not obtain the possible 
driving source of MHO 3284 due to the lack of YSOs in the ambient. 
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Fig. A35.— Same as Fig. Al, but for the region of MHO 3285. Note that we did not obtain the possible 
driving source of MHO 3285 due to the lack of YSOs in the ambient. 
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Fig. A38.— Same as Fig. Al, but for the region of MHO 3288. 


73 









Dec (J2000) 


in n 1 I r 


00 " - . • 


55’ 30" *• 


Q.; • 

.r 


00 " - 


56‘ 30" - 


00 " 


-j-^-r*-1-p-q 


MH03289 




r 


57‘ 30" - 


★ 

37 


-01° 58‘ 00' 


-I_L 


* 1 


_1_L 


18 ^ 32 "^ 14 ® 


10' 

RA (J2000) 

Fig. A39.— Same as Fig. Al, but for the region of MHO 3289. 
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Fig. A40.— Same as Fig. Al, but for the region of MHO 3290. Note that the red circles mark the holes 
inside the polygons. 
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Fig. A41.— Same as Fig. Al, but for the region of MHO 3291. 


76 




Ok (J2000) 



Fig. A42.— Same as Fig. Al, but for the region of MHO 3292. 
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